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ABSTRACT 


a 


H 


The  reduction  of  structural  vibrations  in  ships  and  submarines  is  a  long  standing 
concern  of  the  Navy.  Waveguide  absorbers  are  very  elective  devices  which  can  be  ap¬ 
plied  to  this  problem.  This  study  evaluates  the  increase  in  vibration  damping  of  a  plate 
structure  across  a  broad  frequency  range  using  light  weight  beam  waveguide  absorbers. 
Viscoelastic  and  constrained  layer  beam  waveguide  absorbers  were  studied  both  theore¬ 
tically  and  experimentally.  Impedances  of  the  waveguide  absorbers  at  the  attachment 
point  were  predicted  using  both  Bernoulli-Euler  and  Timoshenko  beam  theory  for  the 
viscoelastic  beam  and  using  sixth  order  beam  theory  for  the  constrained  layer  beam.  The 
theoretically  predicted  impedances  were  compared  with  the  experimental  measurements. 
Results  from  random  vibration  tests  of  a  plate  structure  showed  significant  increases  in 
damping  over  a  broad  frequency  range  (100  Hz  -  2  KHz)  when  the  waveguide  absorbers 
were  attached  on  the  plate. 
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I.  INTRODUCTION 


Suppressing  noise  and  vibrations  of  ships  and  submarines  i<  very  important  tor  the 
Navy  and  has  been  one  of  the  Navy's  long-standing  concerns.  Much  work  has  been  dmie 
on  the  development  of  vibration  control  means  such  as  isolation,  detuning,  viscoelastic 
damping  and  dynamic  absorption.  All  of  these  approaches  have  been  studied  rather 
thoroughly  including  their  range  of  applicability  and  definite  limits  on  their  vibration 
and  noise  reduction  capability  within  acceptable  weight  and  volume  increases. 

Recentlv.  a  waveguide  absorber  concept  has  been  developed  [Ref.  lj.  which  can 
provide  simple  and  highly  effective  vibration  control  over  a  wide  frequency  range.  A 
"Waveguide"  is  a  structure  along  which  vibrational  waves  can  travel,  if  one  end  of  a 
waveguide  is  attached  to  a  vibrating  structure,  some  vibration  energy  will  travel  along 
the  waveguide  in  wave  forms.  If  the  waveguide  is  treated  with  a  high  energy  dissipation 
scheme,  the  damping  of  the  waveguide  causes  the  amplitude  of  the  waves  to  decrease 
as  they  travel  and  the  waveguide  may  be  expected  to  remove  vibrational  energy  from  the 
structure. 

In  previous  studies  by  L'ngar  and  Kurzweil  [Ref.  2]  and  by  Ungar  and  Williams  [Ref. 
3).  semi-infinite  beams  and  exponentially  tapered  semi-infinite  beam  waveguide  absor¬ 
bers  were  studied.  However,  theoretical  prediction  of  the  driving  point  impedances 
showed  wide  discrepancy  from  the  experimental  results  since  the  experiments  were  per¬ 
formed  with  finite  length  beams  and  the  effects  of  viscoelastic  materials  on  beams  were 
not  considered  in  theoretical  calculations. 

In  this  investigation,  two  kinds  of  high  damping  beam  waveguide  absorber  are 
subjected  to  study;  viscoelastic  beam  and  constrained  layer  beam  waveguide  absorbers. 
The  driving  point  impedances  are  theoretically  predicted  for  the  finite  length  beams  of 
highly  viscoelastic  behavior.  Experimental  studies  are  also  performed  on  the  impedances 
of  these  waveguide  specimens  and  on  their  contribution  to  damping  increase  of'  a  test 
plate  structure. 


II.  THEORETICAL  ANALYSIS 

A.  LOSS  FACTOR  CONTRIBUTION  OF  A  WAVEGUIDE  ABSORBER 
Consider  a  linear  structure.  S.  under  a  harmonic  load 

F—f^sinuji  (2.1) 

at  location  c  (Figure  I  on  page  3).  The  structure  5  will  vibrate  with  the  energy  of  vi¬ 
bration 


ir0--L.\/r; 


where  M  represents  the  total  mass  of  the  structure  and  F„  is  the  magnitude  of  average 
velocity.  The  energy  dissipated  in  S  per  cycle.  D0  depends  on  the  original  loss  factor.  ?/„. 
of  the  structure  and  the  following  relation  holds: 

D0 

For  a  lightly  damped  structure.  is  much  smaller  than  1  and  Dn  is  much  smaller  than 
2-lF„.  If  a  point  a  of  a  second  body  B.  a  waveguide  absorber,  is  attached  to  body  5 
at  location  b  of  body  5  .  internal  force,  F will  interact  between  body  S  and  body  B  and 
some  vibration  energy,  D  will  be  dissipated  in  body  B  (Figure  1  on  page  3).  Moreover, 
this  internal  force  inlluer-'es  the  vibrational  motion  of  body  S  and  may  reduce  the  en¬ 
ergy  of  vibration.  IF.  of  body  S  under  the  same  harmonic  load  at  the  same  location  c  . 

I  he  interaction  between  bodies  5  and  B  depends  on  the  impedances  of  bodies  .S  and 
B  at  the  attachment  point.  Z.  and  Z4.  Impedance  is  the  rat’O  of  the  harmonic  force.  F. 
acting  on  a  point  a  of  the  structure  to  the  velocity,  F,  of  the  point. 


and  can  be  describd  by  a  complex  number. 


Z  =  R  +  /.V 


WPG# 


Figure  1.  Vibrating  structure  with  waveguide  absorber. 


L_ 

where  R  is  the  real  part  of  Z  and  resistance  component,  or  simply  resistance,  and  X  is 
the  imaginary  part  of  Z  and  reactive  component,  or  reactance. 

The  energy  loss  per  cycle,  D,  of  a  structure  by  an  attached  waveguide  absorber  is, 

as  shown  by  the  previous  study  [Ref.  2  .  and  Ref.  3], 


where  to  denotes  the  vibration  frequency  in  radian  sec.  V,  represents  the  velocity  of  the 
structure  at  the  attachment  point  of  the  structure  before  the  attachment  of  the  wave¬ 
guide  absorbers,  Z„  and  Z,  denote  the  impedance  of  the  waveguide  absorber  and  the 
structure  at  the  attachment  point,  and  Rb  represents  the  real  part  of  the  impedance  Zy 
From  the  equations  2.2.  2.3  and  2.6.  the  energy  loss  factor.  >?,  due  to  the 
contribution  of  an  attached  waveguide  absorber  can  be  expressed  as  follows. 


Equation  (2.6)  shows  that  the  vibration  energy  loss  ol'a  structure  through  the  u,n- 
eguide  absorber  depends  on  the  magnitude  of  tire  velocity,  E.  at  the  attachment  point 
on  the  structure.  If  this  attachment  point  is  a  node  for  which  V.  =  <>.  then  die  eneig;.  f  o 
factor.  equals  zero.  If  the  attachment  point  is  an  antinode,  then  1  -p—  :  c.m  be  rela¬ 
tively  large.  Equation  (2.7)  indicates  that  the  loss  factor  contribution  i«  smail  il  the  im¬ 
pedance  of  the  beam.  Zk  ,  is  considerably  greater  than  the  impedance  of  the  structure. 
Z,.  If  Z,  is  much  smaller  than  Zs.  then  the  structure  will  not  be  affected  by  the  absorber 
and  the  vibration  energy  of  the  structure,  H'0  in  equation  (2.2).  will  not  be  reduced  much. 
Therefore,  impedance  matching  between  the  structure  and  the  waveguide  absorber  is 
important  to  get  the  highest  energy  loss  factor,  >/.  So,  analytic  methods  which  can  predict 
the  impedance  of  beams  using  the  theory  of  elastic  wave  propagation  in  beams  are  de¬ 
veloped  in  the  next  section. 

B.  IMPEDANCE  OF  BEAM  WAVEGUIDE  ABSORBERS  FROM  WAVE 
PROPAGATION  THEORY 

1.  VISCOELASTIC  BEAM  USING  BERNOULLI -EULER  BEAM  THEORY 

Previous  research  [Ref.  3]  considers  the  impedances  of  the  infinite  viscoelastic 
beam  using  Bernoulli-Euler  beam  theory.  In  this  section,  this  previous  work  is  expanded 
to  include  the  finite  viscoelastic  beam.  The  viscoelastic  beam  will  be  excited  harmonically 
at  its  center  by  a  transverse  force.  F  as  shown  in  Figure  2  on  page  6.  Eor  a  beam  which 
is  free  of  lateral  loading  and  under  an  assumption  that  cross-sectional  areas  remain  plane 
and  normal  to  the  neutral  axis,  the  equation  of  motion  becames  (Ref.  2]. 
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where  yU.t)  is  the  transverse  displacement,  and  E'  is  the  complex  moduli  of  the  beam: 

£*  =  £(  I  +  h/vc).  (2.6) 


and. 

E  :  modulus  of  elasticity  of  the  beam 
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>jLt :  loss  factor  of  the  viscoelastic  material 
m  :  mass  per  unit  length  of  the  beam 
p  :  density  of  the  beam  material 
A  :  cross  sectional  area  of  the  beam 
/  :  moment  of  inertia  of  the  cross  sectional  area 

The  general  form  of  ylx.r)  for  a  propagating  harmonic  wave  in  this  Bernoulli- 
Euler  beam  is 


yU.M  =  Y(>ei{k'x~u,,)  (2.101 

where  A*  is  the  complex  wave  number,  and  oj  is  the  frequency  of  the  propagating  wave. 
Inserting  the  wave  equation  (2.10)  into  the  equation  of  motion  (equation  2.S).  the  dis- 
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persion  relation  is 


(2.11) 


This  equation  gives  the  four 
A*  =  /A(x'l 


ExI(kxf  =  pAu)1 
ditferent  complex  wave  numbers  as  function  of  frequencies. 


tan  Ve  ,  .  .  tan  S/ve  „ 

•)  -  t  sin( - -j - )| 


.2  »-«/•», 


Aq  =  A(N  1  +  >ile )  l,J|ew( 


4 

tan” 

i 

fJ  ve 

4 

| 

tan 

five 

4 

-1 

tan 

*/  ve  . 

_  ^  tan  tan 

<2  =  I  +  five  )  '  (C05( - j - -  )  ~  i  Sin( - ~ 


tan  'tjve 


(2.12.J) 

(2.12.6) 

(2.12.C) 


tan  l>ive 


uui  f I . . j 

—  /  sm( - - — —  )]  (2.12.£f) 


where  subscript  ve  is  an  abbreviation  of  Viscoelastic  beam  using  Bernoulli-Euler  beam 
theory  and  A  is  the  wave  number,  given  by 
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le  cross  sectional  area.  (  ,  — -  )  and  c,  denotes  the 
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r  denotes  the  radius  of  avration  of  the 
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longitudinal  wavespeed  in  the  beam  material  (x  —p  ).  rherelorc,  the  equation  (2.10)  is 

expanded  using  wave  number  k;,  !q.  A;  and  A;. 


y{x.t)  =  ( Yxelk'x  +  Y2elk-'x  +  Y2eik'x  +  Ydeik:x)e~Uot  (2. 14) 


where  Y2,  Ys  and  are  constants  whose  values  depend  on  the  boundary  conditions 
and  .t  is  the  axial  coordinate.  The  complex  transverse  force,  F,  is  denoted  by 


F(.r)  =  £7 
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a.  INFINITE  VISCOELAS  TIC  BE.  IM 

Let's  consider  the  infinite  viscoelastic  beam.  If yfx)  is  to  remain  finite  for  a 
large  x.  then  Yt  must  vanish.  If  there  is  no  wave  coming  toward  the  origin  from  the  po¬ 
sitive  x  direction,  then  K2  must  vanish.  K,  and  Y}  in  equation  (2.14)  are  determined  from 
two  boundary  conditions  at  x  =  0,  as  functions  ofy  and  6,  as  follows: 

>1,0)  =  y  (2.16) 
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At  the  driving  point,  the  input  transverse  force.  /aU).  can  be  caluilatec!  ucing 
the  equation  (2.15).  which  can  be  expressed  using  impedance  and  velocity  at  a=)»  a»> 
follows: 


fl.Ol  =  Z,  r;o)  +  Z2ii(0)  (2.15) 

The  complex  transverse  velocity.  C.  and  the  complex  angular  velocity.  12.  are  related  to 
the  corresponding  displacement,  y  and  0  of  equations  (2.10  and  2.17).  as  follows: 

~  =  "7p*  =  —  iuj  (2.19) 

Z,  represents  the  impedances  of  the  transverse  force  at  the  driving  point  of  the  center 
of  the  infinite  viscoelastic  beam.  If  there  is  no  rotation  at  the  center  of  the  beam,  that 
is  !2(0)  =  0.  the  equation  (2.18)  becomes 

F(0)  =  Zir(U.)  (2.20) 

T  herefore,  the  impedance  of  an  infinite  Bernoulli-Euler  beam  at  the  dining  point.  Z.„. 
can  be  calculated  from 
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where  subscript  ive  is  an  abbreviation  of  the  Infinite  Viscoelastic  beam  using  Betnoul- 
li-Luler  beam  theory. 

b.  FIM  TE  \  ISCOELAS  TIC  BEAM 

For  a  finite  beam,  four  constants  K.  T;  and  }',.  are  determined  from  four 
boundary  conditions  at  x  =  0  and  .v  =  /,  as  follows: 

.iO)  =y  (2.22) 
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From  the  same  considerations  as  those  in  the  previous  section,  the  impe¬ 
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is  determined  as  function  of  frequencies,  where  the  subscript  f:e  is  the  abbreviation  of 
the  Finite  Viscoelastic  beam  using  Bernoulli-Euler  beam  theory. 

2.  FINITE  VISCOELASTIC  BEAM  USING  TIMOSHENKO  BEAM  THEORY 
The  Bernoulli-Euler  beam  theory  was  obtained  using  two  assumptions  which 
neglect  shear  deformation  and  rotational  inertia.  The  Bernoulli-Euler  beam  theory  may 
be  useful  for  low  frequencies  and  long  wavelengths.  However,  the  Timoshenko  beam 
theory  gives  more  accurate  simulation  at  wider  frequency  ranges  and  wave  lengths. 
Therefore  we  will  analyze  the  impedances  of  the  finite  viscoelastic  beam  using  Timosh¬ 
enko  beam  theory  in  this  section.  Two  independent  variables,  the  transverse  displace¬ 
ment.  r(.v.r).  and  angular  displacement.  <}>(xj).  are  defined  in  Timoshenko  beam  theory. 
Beam  motions  are  governed  by  the  following  equations  [Ref.  4[. 
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where  af  =  -Ky  a\ 


—  and  kl  =  -y-  The  complex  moduli  G'  and  E‘  are  defined  as: 


G*  =  G(l  +  iijv!) 


(2.29) 
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( 2.301 


where  subscript  vt  is  an  abbreviation  of  Viscoelastic  beam  using  Timoshenko  beam  the¬ 
ory,  and 

G  :  shear  modulus  of  the  beam 

E  :  elastic  modulus  of  the  beam 

V.,  :  loss  factor  of  the  viscoelastic  material 


.1  :  cross  sectional  area  of  the  beam 
I  :  moment  of  the  inertia  of  the  cross  sectional  area 


:  density  of  the  beam  material 
k  :  shear  deformation  factor 

The  general  forms  yl.v./)  and  dq.v.r)  for  the  propagating  harmonic  \va\e  in  a  Ti¬ 
moshenko  beam  are 

.  .  *  -  Uk  x—ioi)  .  **  i  , 

0{x.i)  —  <P0t 

where  k'  is  the  complex  wave  number,  and  to  is  the  frequency  of  the  propagating  wave. 
Inserting  the  wave  equations  (2.31  and  2.32)  into  the  governing  beam  equations  (2.27 
and  2.2S)  becomes 

(to2  -  a2{k*)2)  T0  +  (  -w,V)<I>0  =  0  (2.33) 

( in |" A' 2 k  )  1q  +  (to"  —  tq  k$  —  u2(^  )~)df)  =  0  (2.c4) 

The  dispersion  relation  of  waves  is  obtained  from  equations  (2.33  and  2.34) 

cq  a2(/c*)J  —  (a2  +  ci2)a>“(Ax)2  +  (to-1  —  c7f  A;,,  aT)  =  0  (2.3o) 


with  the  relationship  between  T0  and  <£q. 
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T  he  equation  (  2.35)  gives  the  four  dilferent  wave  numbers  as  function  oi  wave  frequen¬ 
cies. 
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1  herefore.  the  equations  (2.31  and  2.32)  arc  expanded  using  A,\  A:\  A;  and  Aq 
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The  equation  (2.4ot  becomes  the  following  equation 


<Di.n  =  R{  )\eik'x  +  R2  Y2eik*x  +  R3  Y3c,k*x  +  Rj/klx 
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where 


paG  —  k('i*(  A*)" 
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for  a;  =  1.2.  3.  4. 

For  a  finite  viscoelastic  beam  with  a  length  /.  which  is  free  at  x  =  l  and  excited 
at  .v  =  0,  four  boundary  conditions  are  defined  as: 


v(0)  —  y 


(2.43) 


(MO)  =  0 


(2.44) 


_  c<htl\ 

mu)  =  r:i— —  —  o 

ex 


( 2.45) 


SU)  -  kG  It  -r - </»)  =  0 

<  x 


f  rom  those  boundary  conditions,  the  constants  Y2.  }3  and  Y\  are  determined 

as  functions  of  y  and  0.  At  the  driving  point,  x  =  0.  input  transverse  force. 


/■"(i*)  =  —  5(<M.  can  be  calculated  using  equation  (2. -lot.  which  can  be  expressed  using 
impedances  and  velocity  at  the  driving  point  as  follows; 

/•if)  l  =  Z,  1(0)  +  Z:cb(0|  1 2 .4“  i 


where  the  complex  transverse  velocity.  V,  and  the  complex  angular  velocity.  T>.  are  re¬ 
lated  to  the  corresponding  displacement  as  follows: 

V  <I> 

—  =  —  =  —  i(0  (2.4a) 

>  0 

If  there  is  no  angular  displacement.  0.  at  the  center  of  the  beam,  that  is.  =  u.  the 
equation  (2.47)  becomes 

F(0)  =  Z,  1(0)  (2.4h) 


Therefore,  the  impedance  of  the  Timoshenko  beam  at  the  driving  point.  Z.w,  can  be  cal¬ 
culated  using  the  following  equation. 
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where  subscript  Jvt  is  an  abbreviation  of  the  Finite  Viscoelastic  beam  using  Timoshenko 
beam  theory. 

3.  FINITE  CONSTRAINED  LAYER  BEAM 

In  this  section,  impedances  of  the  constrained  layer  beam,  which  has  viscoelastic 
material  between  the  two  elastic  layers,  are  studied.  The  theory  of  a  damped  sandwich 
beam  was  developed  by  R.A.  Ditaranto  [Ref.  5|.  who  extended  Kerwin's  [Ref.  6  j  bas¬ 
ically  similar  analysis.  D.J.  Mead  and  S.  Markus  [Ref.  7)  expanded  the  concept  by  con¬ 
sidering  these  earlier  works.  The  result  of  their  efforts  is  a  sixth  order  differential 
equation  of  motion  which  is  expressed  in  terms  of'  the  transverse  displacement,  y.  and 
longitudinal  displacement,  u.  for  the  constrained  layer  beam.  The  constrained  laser 
beam  will  be  excited  harmonically  at  its  center  by  a  shear  force.  5.  as  shown  in 
I  igurc  3  on  page  12.  This  equation  of  motion  leads  to  the  impedances  of  the  sandwich 
beams  with  six  boundry  conditions,  which  are  used  for  evaluating  impedances  depending 
on  frequency,  viscoelastic  material  properity,  and  two  elastic  material  properties. 

From  the  |  Ref  7|.  the  governing  sixth  order  differentia!  equation  is 
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Figure  3.  Coordinates  and  sign  conventions  for  a  senn-finite  constrained  layer  beam. 
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where 
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with  the  relationship  between  y(.v.r)  and  u{x.[)  expressed  as: 
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The  parameters  in  this  equation  are  defined  as  follows: 
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where  subscript  c  is  an  abbreviatioti  of  constrained  layer  beam,  and 
G  :  shear  modulus  of  the  viscoelastic  layer 
E  :  elastic  modulus  of  the  elastic  layer 
>lc  :  loss  factor  of  the  viscoelastic  layer 
/  :  moment  of  the  inertia  of  the  cross  sectional  area 
Ii  :  thickness  of  the  constrained  layer  beam 
m  :  mass  per  unit  length  of  the  whole  three  layer  section. 

The  general  form  of  transverse  displacement,  jl.v./h  and  longitudinal  displace¬ 
ment.  for  the  propagating  harmonic  waves  in  a  constrained  layer  beam  are 


y(x.i)  =  T0<? 


i(k\x— col) 


u[x.r)  =  L,')e  ( — '  ) 

where  kr  is  a  complex  wave  number  and  to  is  the  frequency  of  the  propagating  wave. 
After  inserting  the  wave  equations  (2.56  and  2.57)  into  the  sixth  order  beam  equation 
<2.51 1,  the  dispersion  relation  is  obtained  as  follows: 


(/:  Y  +  B,i  I  +  Z)(/<  V  -  r—  (/Yr  -  —77—  =  0 


There  are  six  diiferent  roots  with  six  different  complex  wave  numbers  m  the  above 
equation,  which  means  a  wave  of  six  different  modes  propagate  through  the  cmwi  amed 
layer  beam.  From  the  equations  id. bn  and  2.5”*.  w e  may  rea nance  equation  •' 2.55 » 
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The  general  form  of  the  transverse  displacement  equation  (2.5oi.  y.  and  longi¬ 
tudinal  displacement  equation  (2.5~i.  it  ,  are  obtained  using  the  six  complex  wave  num¬ 
bers  and  a  ratio  Rv  The  equations  [2.5b.  2.57  and  2.59)  as  follow: 


=  ( }\e*'x  +  Y2e-  T  +  Y3e:r>x  +  Y-e  l~  +  YselA,x  +  K-b1'’ dr'1"'  (2,'m  | 

MI.V./I  =  (R{  Yxelk'x  +  R,  Y2e;k:-X  +  R,  Y3e;k*' +  /?a  T4<f ';  v 


1 2.92) 


where 


iBZZDX 

\ k„  I  +Bt.)EiII}Cc 


I  2. he  1 


n  —  1.  2.  3,  4,  5.  6  For  a  constrained  layer  beam  with  a  finite  length  /  which  is  free 
at  .r  -  /  and  excited  at  one  end  at  ,v  =  <>.  the  six  boundary  conditions  are  defined  as  fol¬ 
lows  [Ref.  7|: 
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The  six  constants  Y  (Z  =  l.  2.  3.  4.  5  and  6i  are  determined  from  these  six  boun¬ 
dary  conditions  and  are  a  function  olW,  0  and  u.  At  the  driving  point,  .v  = 66  the  input 
transverse  force.  /3<v)  =  —  Si'M.  moment.  .\Z(<M.  and  normal  foice.  .Ypn.  can  be  calculated 
using  equations  (2.67,  2.6S  and  2.64).  These  input  force  and  moment  can  be  expressed 
using  impedances  and  velocities  and  angular  velocity  at  the  driving  point  as  follows: 


Ft O)  =  Z,  HO)  -j-  ZA2(0)  4-  Z,Zf (0) 


(2.  6] 


Where  the  complex  transverse  velocity.  and  the  complex  angular  velocity. 
£2.  and  the  complex  longitudinal  velocity.  V.  are  related  to  the  corresponding  displace¬ 
ment  as  follows: 


i_  ^  IL  =  I_  =  _  Uo 


(2.71  j 


If  there  is  no  rotational  displacement.  0,  and  no  longitudinal  displacement,  u.  at  the 
center  of  the  beam.  Zy  represents  the  impedances  of  the  transverse  force  at  the  driving 
point  of  the  constrained  layer  beam.  That  is.  £2(0)  =  ()  and  (.’<<>)  =  0.  Therefore,  the 
equation  (2.70)  becomes 


/••(fii  =  z,  noi 


(2.  21 


Therefore,  the  impedances  of  the  constrained  ’aver  beam  at  the  driving  point.  Z .  can 
be  calculated  using  the  following  equation. 
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Where  subscript  c  is  an  abbreviation  ol  the  constrained  la} or  beam. 

C.  COM  PITER  IMPLEMENTATION 

1  he  computer  programs  using  Reduce  and  FORTRAN  languages  are  wed  for  cal¬ 
culation  of  the  theoretical  impedances  from  the  beam  theories  in  this  paper.  Reduce  was 
developed  by  Anthony  C.  Hearn  [Ref.  S),  and  is  a  programming  language  which  solves 
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algebraic  operations  non-numerically.  It  can  manipulate  polynomials  in  a  '.ante  of  fonm 
(substitution,  differentiation  and  integration)  and  solves  one  or  more  simultaneous  al¬ 
gebraic  equations.  After  the  impedances  were  evaluated  bv  a  Reduce  program  uuh 
constant  variables,  a  FORTRAN  program  calculated  the  impedance  with  real  dimen¬ 
sions  and  properties  of  the  viscoelastic  and  constrained  layer  beam  waveguide  absorber 
for  the  0-21)01)11/  frequency  range.  All  constant  variables  of  the  Reduce  program  were 
defined  in  each  FORTRAN  program.  We  tried  to  use  the  same  notations  for  the  Reduce 
program,  the  FORTRAN  program  and  the  text  i  Appendix  A). 

Three  Reduce  programs  and  a  FORTRAN  program  were  used  to  study  each  theory. 
Five  first  Reduce  program  defined  the  moment,  force  and  displacements  which  included 
the  constants  Y...  where  n  equals  1  and  3  for  an  infinite  viscoelastic  beam:  n  equals  l,  2, 
3  and  4  for  a  finite  viscoelastic  beam  and  »  equals  1.  2.  3.  4.  5.  and  6  for  a  constrained 
layer  beam.  F,  can  be  solved  with  boundary  conditions.  After  running  the  first  Reduce 
program,  new  constant  terms  are  defined  for  The  second  Reduce  program  includes 
those  new  constant  variables  which  were  the  results  of  the  first  Reduce  program.  The 


second  Reduce  program  gives  the  constant  variables  finally,  the  third  R.cducc  pio- 
gram  can  evaluate  the  impedance  of  the  beam  with  longitudinal  and  angular  displace¬ 
ment  (and  normal  displacement  for  the  constrained  layer  beam).  The  process  of  the 
constrained  layer  beam  is  similar  up  to  the  second  Reduce  program.  However,  the 
computer  finds  it  hard  to  handle  6  simultaneous  equations  with  6  non-numerical  vati- 
ables,  F„.  Therefore,  we  solved  F4,  Fs  and  F„  in  the  second  Reduce  program  and  then 
substitute  these  values  to  F,.  F,  and  F3  in  the  third  Reduce  program. 

Wave  numbers  as  functions  of  frequencies,  i.e.  dispersion  relation,  of  waves  in  wav¬ 
eguide  absorbers  were  described  in  previous  sections.  It  was  easy  to  find  the  wave  num¬ 
ber  of  a  viscoelastic  beam  through  Bcrnoulii-Fuler  beam  theory  and  Timoshenko  beam 
theory  in  equations  (2.12  and  2.38).  However,  wave  numbers  of  the  constrained  layer 
beam  were  expressed  with  a  complex  six  order  polynomial.  Therefore,  they  were  calcu¬ 
lated  using  Newton’s  method.  The  equation  (2.38)  is  rearranged  using  r*\il  constants 
h.,  r.  and  </,  and  imaginary  constants  h,  and  d.. 


(A.  )  4-  ( br  4-  ib;}{  R  )  -4  <;r(  R  )  4-  ( <ir  4-  ui ;)  =  *• 


(2.M) 


u  here 


A”,  =  (A> 


(2.  /t*.cf) 


h;  =  Inn  IK  4-  Bf.Z) 


—nun 


y 

dr  =  Re{ - jy - ) 


d-t  =  fm( - - - 


i:.' 


12.75/1 


Let  the  complex  variable  K~  =  t  +  iq  (t  and  q  are  real  numbers).  The  equation  (2. ”4) 
becomes  two  equations  coupled  with  two  variables  (t  and  q). 

Jlt.q )  =  t2  +  brt2  +  (  -}q2  -  2 btq  4-  cr)t  +  (  -brq2  +  dr)  ( 2.  "6) 

g(r.q)  =  q 3  +  b-tq2  +  (  -3/2  -  2 brt  -  c,)q  +  {  - b it2  -  dj)  (.2.77 ) 

Starting  from  initial  guesses  (r0,  <7,,)  which  are  near  the  roots  iteration  was  continued  until 
either  the  successive  t  and  q  values  converged  to  the  certain  values  or  the  value  of  the 
function  is  sufficiently  near  zero.  After  a  few  iterations  of  equations  (2.76  and  2.77) 
using  a  FORTRAN  program  (Appendix  A),  we  can  get  three  different  square  of  wave 
number  values.  A'\  were  obtained.  So,  finally,  this  method  calculates  six  different  wave 
numbers.  A;,  from  ±  s  A''  .  where  n  equals  1.  2.  3.  4.  5.  and  6. 

Impedances  of  a  20"  viscoelastic  beam  were  calculated  using  these  computer  pro¬ 
grams  for  viscoelastic  beam  properties  of  200.000  psi  constant  shear  modulus  and  0.2 
energy  loss  factor.  The  results  of  these  simulations  are  shown  in  Figure  4  on  page  19 
and  Figure  5  on  page  20  which  reflect  sharp  resonance  frequencies.  These  sharp  reso¬ 
nance  phenomena  reduced  for  the  same  beam  with  higher  energy  loss  factor  of  0.5.  For 
cases  using  the  Bernoulli-Fuier  beam  theory,  the  impedances  of  the  finite  viscoelastic 
beam  approached  to  the  infinite  viscoelastic  beam  waveguide  absorber  with  increasing 
energy  loss  factor  and  with  increasing  frequency.  The  impedance  using  Timoshenko 
beam  theory  differs  from  the  impedance  from  the  Bernoulli-Fuler  beam  theory 
(1  igure  6  on  page  21  and  Figure  7  on  page  22).  because  ’Timoshenko  beam  theory 
considers  shear  deformation  and  rotary  inertia  effect  of  the  viscoelastic  beam. 
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Impedances  of  a  2n"  constrained  layer  beam  with  viscoelastic  layer  properties  of  2<>u  psi 
constant  shear  modulus  and  O.o  and  l.<*  energy  loss  factors  were  calculated.  Figure  S 
on  page  23  and  Figure  9  on  page  24  compare  the  impedances  vs.  frequency  depending 
on  the  energy  loss  factor  of  the  constrained  layer  beam.  The  impedance  of  the  high 
damping  constrained  layer  beam  shows  the  smooth  and  shifted  resonance  frequencies 
compared  with  the  low  damping  constrained  layer  beam. 
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Figure  5.  Imaginary  part  of  the  16"  viscoelastic  beam  waveguide  absorber  theore¬ 
tical  impedances  [>i  =  0.2)  . 
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Figure  7.  Imaginary  part  of  the  16"  viscoelastic  beam  waveguide  absorber  theore¬ 
tical  impedances  hj  =  0.5i  . 
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Imaginary  part  of  the  2U"  constrained  layer  beam  waveguide  absorber 
theoretical  impedances. 


ML  EXPERIMENT 


A.  IMPEDANCES  OF  THE  TEST  PLAT  E 

A  rectangular  aluminum  plate  with  clamped  boundary  condition  was  selected  as  the 
test  structure,  because  it  has  enough  closely  space  modes  in  the  frequency  range 
t  Ii.'Uf  Lj-2Kz).  An  aluminum  block  ( 30"  x  22"  x  2  )  was  carved  inside  to  make  a  24"  x  lb" 
x  h>. 3  125  '  plate  with  clamp  boundary  condition. 

The  impedances  of  the  test  plate  were  measured  from  impact  hammer  tests.  The 
aluminum  plate  was  excited  by  a  PCB  impact  hammer.  Input  impact  force  was 

measured  by  a  force  transducer  at  the  tip  of  the  impact  hammer.  Responses  were  meas¬ 
ured  by  an  Endevco  2250A-l<>  piezoelectric  accelerometer  which  was  attached  under  the 
impact  point  for  locations  1.  2  and  3  (Figure  10  on  page  26  and  Figure  11  on  page 
2~).  Location  1  is  the  center  of  the  plate.  Location  2  is  the  symmetric  point  of  location 
4.  where  a  shaker  will  be  attached  during  the  damping  measurement  tests.  Location  3  is 
a  general  point  which  is  not  on  any  line  of  symmetry.  A  baseband  measurement  for  each 
location  was  taken  from  0  to  2<X>0Hz.  An  average  of  15  data  samples  was  taken  for  each 
measurement  and  analyzed  using  a  HP-3562A  dynamic  signal  analyzer.  Impedance  of 
the  aluminum  plate  at  each  location  was  determined  after  a  predetermined  data 
processing  procedure  using  the  HP-3562A  dynamic  signal  analyzer. 
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Figure  11.  Test  plate  impedance  measurement  configuration 


B.  IMPEDANCES  OF  THE  W.WEDLTDE  ABSORBER 


Tuc  viscoelastic  beams  and  two  constrained  layer  beams  were  selected  i'  the  wav¬ 
eguide  absorber  specimens.  Dimensions  and  physical  properties  are  shown  in  1"  a  Die  !. 
Mechanical  properties  of  a  viscoelastic  material  depend  on  temperature  and  v. brain  n 
liequeucies.  The  loss  Factors  and  shear  nvodulli  of  the  viscoelastic  material  used  m  tine 
specimens  at  room  temperature  (75 -To  are  shown  in  Figure  12  on  page  do  and 
Figure  13  on  page  31.  These  data  were  obtained  from  Design  Data  Sheets  in  a  h..<ok  by 
Nashil.  et.  ul.  [Ref.  9).  Analytic  expressions  of  the  loss  factors  and  shear  moduli  were 
determined  from  curve  Fits  to  these  data.  The  following  four  equations  represent  the  Iocs 
factors  and  shear  moduli  of  the  LD-doo  and  ISD-1 12  viscoelastic  material. 
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f  represents  the  frequency.  ?/„  and  >ic  represent  the  loss  factor  of  LD-400  and  ISD-112 
viscoelastic  material.  G„  and  Gc  represent  the  shear  modulus  of  LD-400  and  ISD-112 
viscoelastic  material.  Figure  14  on  page  32  shows  waveguide  absorbers  which  were  used 
in  this  study. 

Impedances  of  waveguide  absorbers  at  the  attachment  points  were  determined  liom 
steady  state  random  vibration  tests.  The  waveguide  absorber  was  excited  by  a  Wilcoxon 
F4  F'7  vibration  generator  using  a  band  limited  white  noise  signal  (o-20oollz)  from  a 
IIP-3562A  dynamic  signal  analyzer  (  Figure  15  on  page  33).  The  waveguide  absorber  was 
mounted  on  the  Wilcoxon  F4  F7  vibration  generator  (Figure  16  on  page  34).  Force  and 
acceleration  were  measured  using  a  force  transducer  and  a  piezoelectric  accelerometer  in 
the  impedance  head  of  the  Wilcoxon  14  17  shaker.  These  signals  were  analyzed  by  a 
IIP-3562A  dynamic  signal  analyzer,  loo  data  samples  were  averaged  for  each  measure¬ 
ment. 

I  he  measurement  of  impedances  included  the  impedances  of  bolts  and  nuts  and 
connecting  aluminum  pieces.  Therefore,  the  impedances  of  bolts,  nuts,  and  connecting 
aluminum  pieces,  the  so-called  mass  elfeet,  were  measured  separately.  Total  impedances 
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C.  DAMPING  OF  THE  PLATE 


Modal  damping  values  o!"  the  plate  were  determined  from  structural  transfer  func¬ 
tions  at  a  driving  point  and  between  input  and  output  points  on  the  test  plate  which 
were  measured  using  random  excitations.  The  IIP-35o2A  dynamic  signal  anaiy/er  was 
used  in  the  measurement  of  the  transfer  functions.  The  Wileoxon  13  vibration  generator 
excited  the  plate  at  location  4  using  random  oscillation  signals  From  a  HP-35t>2.\  dy¬ 
namic  signal  analyzer.  Input  forces  and  output  responses  at  the  driving  point  were 
measured  using  a  force  transducer  and  a  piezoelectric  accelerometer  in  the  impedances 
head  of  Wileoxon  F3  shaker  (Figure  17  on  page  3M.  The  damping  measurement  con¬ 
figuration  is  shown  in  Figure  IS  on  page  37.  It  was  necessary  to  measure  responses  at 
several  points  other  than  the  driving  point  to  accurately  measure  damping  values  of 
closely  spaced  modes.  These  were  measured  using  the  Endevco  2250A-I0  piezoelectric 
accelerometer.  Modal  damping  values  of  all  the  modes  in  the  frequency  range  of 
1 1 )0_ 2* M.)Ul iz  were  determined  from  the  response  functions  using  the  Half  power  Band 
Width  method  [Ref.  1  ( )|  and  a  curvefit  method  [Ref.  1 1] 

Modal  damping  values  of  the  plate  without  a  waveguide  absorber  were  measured 
as  the  baseline.  A  waveguide  absorber  was  attached  to  the  plate  at  location  1.  The 
waveguide  absorber  was  changed  after  each  measurement  of  the  damping  for  four  dif¬ 
ferent  waveguide  absorbers  which  were  described  in  section  4-B.  After  the  damping  was 
measured  at  location  1.  measurements  were  taken  at  location  2  and  location  3. 

The  most  ellective  waveguide  absorber  at  each  location  was  selected.  The  20"  vis¬ 
coelastic  beam  waveguide  absorber,  the  16"  constrained  layer  beam  waveguide  absorber 
and  the  16"  viscoelastic  beam  waveguide  absorber  were  most  effective  at  location  1.  2 
and  3.  respectively.  Therefore,  the  damping  with  two  waveguide  absorbers  attached  at 
locations  1  and  2.  and  the  damping  with  three  waveguide  absorbers  at  locations  1.  2  and 
3  were  also  measure,  to  see  the  effect  of  multiple  waveguide  absorbers. 
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IV.  RESULTS  AND  DISCUSSIONS 


A.  IMPEDANCES  OF  THE  ITS  T  PLATE 

Impedances  of  the  test  plate  were  measured  from  input  hammer  tests  as  described 
in  Section  3. A.  Since  the  aluminum  test  plate  is  a  lightly  damped  structure,  the  impe¬ 
dances  at  resonance  frequencies  are  very  low  and  at  anti-resonance  frequencies  arc  \ery 
high.  Phases  with  lSiT  phase  shift  are  near  ±9ir-’  both  at  the  icsonance  and  a*  the 
anti-resonance  frequencies. 

Since  location  1  is  the  center  point  of  the  plate,  only  small  numbers  of  icsonance 
and  anti-resonance  frequencies  appear  (Figure  IP  on  page  39  and  Figure  2o  on  page 
40).  Neither  location  2  nor  location  3  are  on  any  line  of  symmetry  of  the  test  plate:  all 
the  resonance  and  anti-resonance  points  appear  in  the  impedance  vs.  frequency  dia¬ 
grams.  (Figure  21  on  page  41.  Figure  22  on  page  42.  Figure  23  on  page  43  and 
Figure  24  on  page  44).  However,  the  anti-resonance  frequencies  at  location  2  and  those 
at  location  3  are  dilferent. 
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B.  IMPEDANCES  OF  THE  WAVEGUIDE  ABSORBER 

Two  different  kinds  of  waveguide  absorbers,  a  viscoelastic  beam  and  a  constrained 
laver  beam,  are  investigated  in  this  study.  In  a  previous  study  by  Incur  and  Williams 
j Ref.  3|.  impedances  were  predicted  only  for  the  infinite  viscoelastic  beam  usng  Ber- 
noulii-Fuler  beam  theory  and  the  theoretical  prediction  sb.owed  wide  dill’ercnce  from  the 
experimental  results.  In  this  study,  this  previous  work  is  extended  to  include  the  impe¬ 
dances  for  the  finite  viscoelastic  beam.  Since  Bernoulli-Euler  beam  theory  does  not  in¬ 
clude  the  shear  deformation  and  rotary  inertia  effect  of  the  viscoelastic  beam,  the 
predictions  from  the  Timoshenko  beam  theory  aie  studied  to  investigate  these  effects  on 
wave  propagation  and  impedances,  and  are  compared  with  the  results  lrom  the  Ber¬ 
noulli-Euler  beam  theory.  The  impedances  of  the  constrained  layer  beam  waveguide  ab¬ 
sorber  were  predicted  using  the  sixth  order  beam  theory. 

The  impedances  from  these  theoretical  predictions  are  compared  with  tire  impe¬ 
dances  from  the  experimental  measurements  for  the  16"  and  20"  length  waveguide  ab¬ 
sorbers  of  two  different  kinds.  In  the  calculation  of  theoretical  predictions,  the  loss 
factors  and  shear  moduli  of  viscoelastic  material  were  varied  as  function  of  frequency 
as  described  in  Section  3.B.  During  the  experiment,  the  Wikoxon  F4  F7  vibration  gen¬ 
erator  produced  characteristic  errors  which  occured  at  different  frequencies  depending 
on  the  mass  of  the  waveguide  absorbers.  Therefore  the  viscoelastic  beam  has  unreliable 
data  around  7(X>Hz  and  the  constrained  layer  beam  has  unreliable  data  around  35UI  Lt. 

Figure  25  on  page  47  and  Figure  26  on  page  4S  show  the  three  theoretical  impe¬ 
dances  compared  with  experimental  data.  The  Bernoulli-Euler  beam  theory  for  finite 
beams  shows  some  standing  wave  effect  at  low  frequency  range  due  to  reflection  at  the 
free  end  of’ the  beam,  but  approaches  the  results  of  the  infinite  beam  as  frequency  in¬ 
creases.  T  he  results  of  the  Bernoulli-Euler  beam  theory  do  not  quite  correspond  to  ex¬ 
perimental  results  for  the  viscoelastic  beam.  Impedances  from  the  Timoshenko  beam 
theory  closely  follow  the  experimental  results  for  viscoelastic  beams  as  shown  in 
Figure  25  on  page  4~  and  Figure  2(>  on  page  48.  These  show  that  the  shear  deformation 
and  rotary  inertia  effects  are  significant  for  the  impedance  prediction  of  beam  type 
waveguide  absorbers. 

Figure  2~  on  page  44  and  Figure  28  on  page  50  represent  the  real  and  imaginary 
parts  of  the  In"  constrained  la;  -r  beam  waveguide  absorber  and  Figure  24  on  page  51 
and  Figure  3d  on  page  52  represent  the  real  and  imaginary  parts  of  the  2d"  constrained 
laver  beam  waveguide  absorber.  In  both  beams,  the  theoretical  results  of  the  impedances 
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using  the  sixth  order  beam  theory  and  the  resonance  frequency  patterns  are  almost  the 
same  as  the  experimental  result. 

Both  Timoshenko  beam  theory  and  the  sixth  order  beam  theorv  predict  resonance 
frequencies  shifted  amount  from  the  experimental  resuirs  as  frequency  increases,  and  the 
impedance  magnitudes  of  experimental  results  are  bigger  than  the  theoretical  results. 
These  difference  between  theoretical  and  experimental  impedances  may  derive  from  the 
inaccuracy  of  shear  modulus  and  loss  factors  which  values  were  selected  from  75'  F  room 
temperature  and  changed  depending  on  frequency.  However,  the  theoretical  prediction 
using  l  imoshenko  beam  theory  for  the  viscoelastic  beam  and  the  sixth  order  beam  the¬ 
ory  for  the  constrained  layer  beam  follow  the  same  trend  1 1  igure  12  on  page  3<>  and 
Figure  1 3  on  page  3 1 ). 

Impedances  of  wavign'de  absorbers  are  very  low  and  do  not  show  any  sharp  vari¬ 
ation  for  all  frequencies  in  the  range  tested,  which  are  quite  different  from  the  impedance 
of  the  test  plate  of  very  lightly  damped  structures. 
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C.  DAMPING  OF  THE  PLATE 


The  purpose  of  this  experiment  is  to  investigate  the  effect  of  light  weight  '.va\eeuide 
absorbers  on  damping  increase  in  plate  t> pe  structures  over  a  wide  frequency  ranee. 
Damping  of  the  test  plate  structure  was  measured  with  an  attachment  of  each  one  of 
four  different  waveguide  absorbers  (16"  or  20"  viscoelastic  or  constrained  laser  beam i 
at  three  different  locations.  The  frequency  response  and  the  modal  damping  values  of  the 
aiuminum  test  plate  with  a  waveguide  absorber  were  compared  with  those  of  the  same 
plate  without  any  waveguide  absorber  over  a  wide  frequency  range  (O-dOooH/f 

The  modal  damping  values  of  the  test  plate  without  a  waveguide  absorber  are  shown 
in  Figure  31  on  page  56  and  Figure  32  on  page  57.  They  show  that  more  than  2o  ligiitiv 
damped  modes  are  in  the  frequency  range  between  1  <  *< '  1 1  z  to  2,J(.l<*IIz  with  the  lowest 
frequency  1451 \z. 

The  damping  contribution  of  the  same  waveguide  absorber  to  a  same  plate  was  dif¬ 
ferent  depending  on  the  location  of  the  waveguide  absorber  attachment  point.  For  the 
20"  viscoelastic  beam  waveguide  absorber,  there  was  marked  damping  increase  for 
modes  lower  than  1  looHz  when  it  was  attached  at  location  1  (Figure  33  on  page  58  and 
Figure  34  on  page  54).  However,  it  was  more  effective  for  modes  higher  than  1200llz 
when  it  was  attached  at  location  2  (Figure  35  on  page  60  and  Figure  36  on  page  61). 
When  it  was  attached  at  location  3  its  contribution  to  plate  damping  increase  was  ef¬ 
fective  for  wider  frequency  range  of  250Hz-l500Hz  (Figure  37  on  page  62  and 
Figure  3S  on  page  63).  This  trend  applies  to  other  waveguide  absorbers,  though  the 
magnitude  of  damping  contribution  and  the  frequency  ranges  were  dilferent  for  each 
case  i  Figures  in  Appendix  13). 

Different  modal  damping  effects  with  a  dilferent  waveguide  absorber  at  the  same 
location  were  also  observed.  At  location  1.  the  20"  viscoelastic  beam  waveguide  absorber 
was  most  effective  for  the  frequency  range  lower  than  1 100Hz  (Figure  33  on  page  58 
and  1  igtue  34  on  page  54)  and  the  16"  constrained  layer  beam  waveguide  absorber  was 
most  cfTccti v e  for  the  frequency  range  higher  than  I20uHz  (Figure  34  on  page  64  and 
I  igurc  4n  on  page  65).  At  location  2.  the  20"  and  16"  viscoelastic  beam  waveguide  ab¬ 
sorbers  showed  moderate  contribution  to  plate  damping  and  the  2o  constrained  layer 
beam  waveguide  absorber  contributed  the  least  amount  to  plate  damping.  For  frequency 
ranges  over  12'Mill/.  the  16  constrained  layer  beam  waveguide  absorber  produced  the 
largest  damping  increase.  At  location  3.  the  20  and  1 6"  viscoelastic  beam  waveguide 
absorbers  and  the  20"  constrained  layer  beam  waveguide  absorber  produced  moderate 


d  imping  increase  over  a  wide  frequency  range  of  25oH/-l5oi'»I  \i.  However,  the  In 
constrained  laver  beam  waveguide  absorber  showed  the  least  damping  increase. 

Alter  the  selection  of  optimum  damping  at  each  location,  two  wane  guide  absorbers 
>2"  viscoelastic  and  lo"  constrained  layer  beam)  were  attached  at  locations  I  and  2  and 
three  waveguide  absorbers  ( 2o"  viscoelastic,  lo"  constrained  layer  and  !(>'  vhcoekn-tic 
beam)  were  attached  at  locations  1,  2  and  3  to  see  the  effect  of  multiple  w  aveguide  ab¬ 
sorbers  on  the  test  plate  damping  (Figure  41  on  page  6o>.  Figure  42  on  page  07  and 
Figure  43  on  page  6S  represent  frequency  response  and  modal  damping  of  the  plate  with 
two  waveguide  absorbers.  Figure  44  on  page  69  and  Figure  45  on  page  TO  represent 
frequency  response  and  modal  damping  of  the  plate  with  three  waveguide  absorbers. 
Damping  of  the  plate  increased  with  increasing  number  of  the  waveguide  absorbers 
(  Figure  42  on  page  67  and  Figure  44  on  page  69).  The  damping  increase  resulting  from 
attachment  of  two  or  more  absorbers  resemble  the  sum  of  the  seperate  result  of  each. 

Impedance  of  the  test  plate  at  location  1  is  compared  with  the  driving  point  impe¬ 
dance  of  the  2<)"  viscoelastic  beam  in  Figure  46  on  page  71.  This  figure  shows  that  at 
resonance  frequencies  the  impedance  of  the  test  plate  is  much  smaller  than  that  of  the 
waveguide  absorber.  Figure  47  on  page  72  and  Figure  4S  on  page  73  show  the  com- 
parsion  between  the  impedances  of  the  test  plate  at  location  2  and  location  3  with  the 
drving  point  of  the  16"  constrained  layer  beam  and  16"  viscoelastic  beam,  respectively. 
These  figures  show  that  at  low  frequency  range  (<40011?’)  the  impedances  of  the  test 
plate  were  higher  than  those  of  waveguide  absorbers  and  at  high  frequency  range 
(  >40(.)Hz)  impedances  of  waveguide  absorbers  were  higher  at  some  resonance  frequen¬ 
cies. 

from  equation  (2.7),  the  loss  factor  contribution,  >/,  of  a  waveguide  absorber  is 
proportional  to  the  following  3  terms: 
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The  first  term  is  determined  from  waveguide  characteristics  only  and  the  second  term  is 
determined  from  the  impedance  ratio  between  the  test  structure  and  a  waveguide  ab¬ 
sorber  at  the  attachment  point.  However,  the  third  term  depends  on  the  interaction  be¬ 
tween  the  test  structure  and  a  waveguide  absorber  and  the  whole  response  of  the  test 
plate.  If  Z,  is  much  smaller  than  Z.  then  H_  will  remain  constant  before  and  after  the 
attachment  of  a  waveguide  absorber.  Previous  studies  [Ref.  2  and  Ref.  3)  were  based 
on  this  assumption  and  gave  misleading  indications  that  experimental  loss  factors  would 
be  much  higher  than  theoretical  predictions.  In  the  present  study,  as  shown  in 
Figure  4o  on  page  ”1.  Figure  4“  on  page  ”2  and  Figure  4S  on  page  ~3.  the  impedances 
of  the  test  plate  were  similar  in  magnitude  to  impedances  of  waveguide  absorbers  at  re¬ 
sonance  frequencies  and  the  loss  factor  contribution  due  to  waveguide  absorbers  could 
not  be  predicted  using  equation  ( 2. ~  >. 
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MODAL  DAMPING  VS.  FREQUENCY 


Figure  3 


Modal  damping  factors  vs.  frequency  of  the  test  plate  with  a  20"  vis¬ 
coelastic  beam  waveguide  absorber  at  location  1  and  without. 
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MODAL  DAMPING  VS.  FREQUENCY 


Figure  39.  The  driving  point  frequency  response  of  the  test  plate  with  a  16  con¬ 
strained  layer  beam  waveguide  absorber  (dashed)  at  location  1  and 
without  (solid). 


Figure  41. 


The  test  plate  with  three  waveguide  absorbers 


Figure  42.  The  driving  point  frequency  response  of  the  test  plate  with  a  20"  viscoe¬ 
lastic  (at  location  1)  and  16"  constrained  layer  (at  location  2)  beam 
waveguide  absorber  (dashed)  and  without  (solid). 


MODAL  DAMPING  VS.  FREQUENCY 


Figure  43.  Modal  damping  factors  vs.  frequency  of  the  test  piate  with  a  2D"  viscoe 


lastic  (location  l)  and  a  16"  constrained  layer  (location  2)  beam  wave 


guide  absorber  and  without. 
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MODAL  DAMPING  VS.  FREQUENCY 


Figure  45.  Modal  damping  factors  vs.  frequency  of  the  test  plate  v.ith  a  2D”  viscoe¬ 
lastic  (location  1).  a  16"  constrained  layer  (location  2)  and  a  16''  viscoe- 
iatic  beam  waveguide  absorber  and  without. 


V.  CONCLUSIONS 


r 


'['he  purpose  of"  this  study  was  to  develop  two  kinds  of  high  damping  beam  wave¬ 
guide  absorbers  for  application  to  naval  structures  o\er  a  wide  frequency  range 
( l00Hz-2<>U0Hz)  and  to  evaluate  theoretical  prediction  schemes  using  various  beam 
theories.  Comparison  between  theoretical  predictions  and  experimental  results  show  that 
Iimoshenko  beam  theory  and  the  six  order  beam  theory  can  be  used  in  the  prediction 
of  impedances  of  viscoelastic  beam  waveguide  absorbers  and  constrained  layer  beam 
waveguide  absorbers,  respectively. 

Application  of  waveguide  absorbers  to  a  test  plate  structure  showed  that  damping 
of  the  test  structure  can  be  increased  significantly  with  one  waveguide  absorber  over  a 
frequency  range.  The  magnitude  and  the  frequency  range  of  damping  increase  depend 
on  the  impedance  of  the  waveguide  absorber  and  on  the  location.  The  prediction  of 
damping  increase  due  to  an  attached  waveguide  absorber,  equation  (2.7).  does  not  show 
the  effect  of  waveguide  absorber  completely  since  this  equation  does  not  include  the 
structure  vibration  reduction  effects  due  to  the  waveguide  absorber.  Since  the  contrib¬ 
ution  of  each  waveguide  absorber  to  the  damping  of  the  test  plate  is  additive,  the  modal 
damping  values  of  all  the  vibration  modes  of  the  test  plate  in  a  wide  frequency  range  can 
be  increased  using  multiple  waveguide  absorbers. 
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VJ.  RECOMMENDATIONS 

From  the  results  of  present  study  following  are  recommended  for  future  studies  on 
waveguide  absorbers: 

1.  Investigation  of  the  effect  of  a  waveguide  absorber  on  the  reduction  in  vibration 
of  the  original  structure  and  derivation  of  the  relation  between  the  impedances  and  the 
mode  shape  of  the  structure. 

2.  Theoretical  studies  on  impedance  prediction  of  waveguides  of  different  shape  such 
as  circular  viscoelastic  plates  and  application  to  damping  increase. 

3.  Investigation  of  interaction  between  waveguide  absorber  for  multiple  waveguide 
absorbers  applications,  and 

-f  Optimization  of  application  of  waveguide  absorbers  in  design,  number,  and 
location. 
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APPENDIX  A. 


1.  THE  IMPEDANCE  OF  THE  INFINITE  VISCOELASTIC  BEAM  USING 
BERNOULLI- El'LER  BEAM  THEORY 


THE  THREE  REDUCE  PROGRAMS  WERE  USED  FOR  EVALUATING  THE  IMPEDANCE  * 
OF  THE  INFINITE  VISCOELASTIC  BEAM  USING  BERNOULLI -EULER  BEAM 
THEORY. 


C 

C  THE  FIRST  REDUCE  PROGRAM  DEFINED  TRANSVERSE  FORCE  AND  TWO  BOUNDARY 
C  CONDITIONS. 

C 

YY:  =Y1*EXP( I*K1*X)+Y3*EXP( I*K3*X); 

PP:  =DF ( YY , X ) ; 

MM:  =-EI*DF(PP,X); 

FF:  =-DF(MM,X); 

YO:  =SUB(X=0,YY)-Y; 

PO: =SUB(X=0,PP)-P; 

BYE; 

C 

C  THE  SECOND  REDUCE  PROGRAM  SOLVED  CONSTANT  VARIABLES  (Y1  AND  Y3). 

C 

Y0:=Y1+Y3-Y; 

P0:=P1*Y1+P3*Y3-P; 

SOLVE ( LSTC YO , PO) ,Y1,Y3); 

BYE; 

C 

C  THE  THIRD  REDUCE  PROGRAM  EVALUATED  THE  TRANSVERSE  FORCE. 

C 

YY:  =Y1*EXP( I*K1*X)+Y3*EXP( I*K3*X); 

FFF:  =E I *DF ( YY , X , 3 ) ; 

FF: =SUB(X=0 ,FFF); 

BYE; 


,  jl  j.  . 
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THIS  FORTRAN  PROGRAM  DEVELOPED  THE  IMPEDANCE  OF  THE  INFINITE 
VISCOELASTIC  BEAM  USING  BERNOULLI -EULER  BEAM  THEORY. 

DEFINED  VARIABLES; 

COMPLEX  YOUNG'S  MODULUS  OF  THE  VISCOELASTIC  BEAM 
FREQUENCY  (HZ) 

TRANSVERSE  FORCE  ( LBF) 

COMPLEX  SHEAR  MODULUS  OF  THE  VISCOELASTIC  BEAM 
GRAVITY  (  32316  IN/SEC 
IMPEDANCE  OF  THE  VISCOELASTIC  BEAM 
POISSON'S  RATIO 

DENSITY  OF  THE  VISCOELASTIC  BEAM 


E 

F 

FF 

G 

GRAV 

IMP 

POISS 

ROH 


XB,  XH , XL  :  DIMENSION'S  OF  THE  BEAM  (WIDTH,  HEIGHT,  LENGTH) 
XI  :  AREA  MOMENT  OF  INERTIA 

XK1 ,  XK2,  XX 3 ,  XK4  :  WAVE  NUMBERS 

XNETA  :  ENERGY  LOSS  FACTOR  OF  THE  VISCOELASTIC  BEAM 
VI  :  IMAGINARY  (0,1) 


C 

COMPLEX  VI ,G ,E ,KB 1 ,KB2 ,XK1 ,XK2 ,XK3 ,XK4 , PI ,P3 , A1 , A2 , A3 , A4 ,FF , IMP 
REAL  F , POI SS , GRAV , P I , XB , XH , XG , ROH , XI 
REAL  W , XLOG , XNETA , XA , MCA , XSA , KC , IMPR , IMP I 
C 

OPEN(UNIT=15 ,FIL£=' IFI ' ) 

C 

C  CONSTANTS 
C 

F=10. 0 
POISS=0. 45 
GRAV=386. 0 
PI=3. 1415926 
V1=CMPLX( 0.  0,1.  0) 

C 

C  DIMENSIONS 

C 

XB=1.  0 
XH=0.  35 

L» 

C  PROPERITIES 

C 

ROH=0.  05505/GRAV 
XI=( XB*XH**3 ) / 12.  0 
C 

C  FOUR  DIFFERENT  WAVE  NUMBERS 

C 

DO  100  1=1,100 
W=2*PI*F 

XL0G=ABS( L0G10( F/62 ) ) 

XNETA=0.  65*EXP( -0.  52732*XLOG**l.  956) 

C  XNETA=0.  2 

C  XNETA=0.  5 

XG=0. 0000250  3*F** 3 - 0 . 1752*F**2+457. 5883*F+29280 
C  XG=200000 

G=XG*CMPLX( 1.  0 , XNETA) 

E=2*G*( 1+POISS) 

XA=ATAN( XNETA) / 4. 0 
XCA=COS( -XA) 

XSA=SIN( -XA) 

KC=( S0RT( 1+XNETA**2) )**( -0.  25 ) 

KB 1=-CSQRT( ROH*XB*XH*W**2/f E*XI ) ) 

KB2=  CSORT(  ROH*XB*XIl*W**2 / (  E*XI ) ) 

XK1=  CSQRTf KB  1 )*KC*CMPLX( XCA.XSA) 

XK2=-CSQRT(  KB  l)*KC*CtJPI.X(  XCA.XSA) 

XK3— -CSQRT( KB  2 ) * K C*C MP LX ( XCA.XSA) 

XK4=  CSQRT( KB 2) *KC*CMPLX( XCA , XSA) 

C 


C  DEFINED  CONSTANTS  FROM  THE  THREE  REDUCE  PROGRAMS 
C 

P1=V1*XK1 
P2=V1*XK2 
F3=V1*XK3 
F4=V1*XK4 
Al=  P3/CP3-P1) 

A2=0 

A3=-P1/(P3-P1) 

A4=0 

FF=- V 1*E*XI*( A 1*XK 1**3+A3*XK3**3 ) 

C 

C  IMPEDANCES 
C 

IMP=2*FF/(V1*W) 

IMPR=REAL( IMP) 

IMP I =A I MAG (IMP) 

WRITE( 15 , 1000)  F , IMPR, IMP I 
F=F+20.  0 
100  CONTINUE 

1000  FORMAT ( IX, D9. 3,2X,D12. 5,2X,D12. 5) 

CL0S£(UNIT=15) 

STOP 

END 

2.  THE  IMPEDANCE  OF  THE  FINITE  VISCOELASTIC  BEAM  USING 


=Y1*£XP( I *K 1*X ) +Y 2*EXP ( I*K2*X)+Y3*EXP( I*K3*X)+Y4*EXF( I*K4*X) ; 


BERNOULLI-EULER  BEAM  THEORY 

YY 

PP:  =DF ( YY , X ) ; 

=-EI*DF( PP ,X) ; 

=-DF(MM,X); 

=SUB(X=0 , YY) - Y; 

=SUB("  X=0  ,  PP)  -P; 

=SUB(X=L,MM); 

=SUB(X=L,FF); 


MM 

FF 

YO 

FO 

ML 

FL 

BYE 

C 

C 

C 

C 

YO 

PO 

ML 

FL 


THE  SECOND  REDUCE  PROGRAM  SOLVED  CONSTANT  VARIABLES  (Yl,  Y2 ,  Y3 
AND  Y4) . 

=Y 1+Y2+Y3+Y4-Y; 

=P 1*Y 1+P2*Y2+P3*Y3+P4*Y4 - P; 

=M 1 * Y 1 +M2* Y2+M3*  Y 3+M4*Y4 ; 


=F1*Y1+F2*Y2+F3*Y3+F4*Y4; 

SOLVE!  LST(  YO  ,  PO  ,!1L,FL) ,  Yl ,  Y2  ,  Y3  .  Y4); 

BYE; 

C 

C  THE  THIRD  REDUCE  PROGRAM  EVALUATED  THE  TRANSVERSE  FORCE. 

C 

YY:  =Y1*EXP(  I*Kl*X)+Y2*EXFf  I*K2*X)+Y3*EXP(  I*K3*X)+Y4*EXP(  I*K4'AX) ; 
FFF:=EI*DF(YY,X,3); 
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THIS  FORTRAN  FROGRAM  DEVELOPED  THE  IMPEDANCE  OF  THE  FINITE 
VISCOELASTIC  BEAM  USING  BERNOULLI -EULER  BEAM  THEORY. 


DEFINED  VARIABLES; 

E  :  COMPLEX  YOUNG'S  MODULUS  OF  THE  VISCOELASTIC  E 

F  :  FREQUENCY  (HZ) 

FF  ;  TRANSVERSE  FORCE  (LBF) 

G  :  COMPLEX  SHEAR  MODULUS  OF  THE  VISCOELASTIC  BE* 

GRAV  :  GRAVITY  (32)?  IN/SEC 

IMP  :  IMPEDANCE  OF  THE  VISCOELASTIC  BEAM 

POISS  :  POISSON'S  RATIO 

ROII  ;  DENSITY  OF  THE  VISCOELASTIC  BEAM 

XB,  XH.XL  :  DIMENSIONS  OF  THE  BEAM  (WIDTH,  HEIGHT,  LENT 

XI  ;  AREA  MOMENT  OF  INERTIA 

XK1 ,  XK2 ,  XK3 ,  XK4  :  WAVE  NUMBERS 

XNETA  :  ENERGY  LOSS  FACTOR  OF  THE  VISCOELASTIC  BEAM 
VI  :  IMAGINARY  (0,1) 


VpVryrVe^VVcyr*ypVrVpVcVc*VrVrVp*Vc*****VrVfVfiViVVrVrVriVVrVrVryfVr:VVrVrVc,5VVrVrVrVf,}V*VrVr,:V*Vr-,V:VVrVr’;V*-;V'/rtVVr“/c,-5V*-,f-,V 

c 

COMPLEX  VI ,G,E ,KB1 ,KB2 ,XK1 ,XK2 ,XK3 ,XK4 , PI , P2 , P3 , P4 
COMPLEX  Ml ,M2 ,M3 ,M4,F1 , F2 ,F3 ,F4 ,DEN . A1 , A2 , A3 , A4,FF , IMP 
REAL  F, POISS, GRAV, PI, XB , XL, XH,XG,ROH, XI, W,XLOG, XNETA, XA,XCA,XSA 
REAL  KC ,KR1 , KR2 ,KR3 ,KR4 ,KI 1,KI2,KI3,KI4, IMPR, IMP I 
C 

0PEN(UNIT=15 ,FILE=' FFI ' ) 

C 

C  CONSTANTS 
C 

F=10.  0 
P0ISS=G.  45 
GRAV=386.  0 
PI=3. 1415926 
V1=CMPLX(0. 0,1. 0) 

C 

C  DIMENSIONS 

C 

XB=1.  0 
XL=3. 0 
XH=0. 35 
C 

C  PROPERITIES 

C 

R0H=0.  05505/GRAV 
C 

C  FOUR  DIFFERENT  WAVE  NUMBERS 


X I = ( X  B  *XH** 3 ) / 1 2 . 0 
DO  100  1=1,100 
W=2*PI*F 

XLOG=ABS( LOG10( F/62) ) 

XNETA=0. 65*EXP( -0.  52732*XLOG**l.  956) 
XNETA=0. 2 
XNETA=0. 5 


XG=0 .  000025 03*F**3-0. 1752*F**2+457. 5883*F+29280 

C  XG=200000 

G=XG*CMPLX( 1.  O.XNETA) 

E=2*G*( 1+FOISS ) 

XA=AT AM ( XNETA ) / 4 . 0 
XCA=COS( -XA) 

XSA=SIN( -XA) 

KC=( SQRT( 1+XNETA**2 ) )**( -0. 25 ) 

K3 1=-CSQRT( R0H*XB*XH*W**2 / ( E*XI ) ) 

KB2=  CSQRT( R0H*XB*XH*W**2/(E*XI ) ) 

XK1=  CSQRT( KB l)*KC*CMPLXf XCA ,XSA ) 

XK2=-CSQRT( KB  1 )*KC*CMPLX( XCA . XSA ) 

XK3=-CSQRT( KB2)*KC*CMPLX( XCA , XSA ) 

XK4=  CSOST( KB  2 ) *KC*CMPLX( XCA , XSA ) 

KR1=REAL( XK1 ) 

KR2=REAL( XK2 ) 

KR3~REAL(XK3) 

KR4=REAL( XK4) 

KI 1=AIMAG( XK1 ) 

KI2=AIMAG(XK2) 

K 1 3  =A I MAG ( XK 3 ) 

K I 4=A I MAG ( XK4 ) 

C 

C  DEFINED  CONSTANTS  FROM  THE  THREE  REDUCE  PROGRAMS 

C 

P1=V1*XK1 

P2=V1*XK2 

P3=V1*XK3 

P4=V1*XK4 


C 


\Jq 


=EXP( -XL*KI 1 )*(C0S(XL*KR1)+V1*SIN( XL*KR1 ) )*XK1**2 
=EXP( -XL*KI 2 )*( C0S(XL*KR2)+V1*SIN( XL*KR2) )*XK2**2 
=EXP( -XL*K 1 3 ) ' *( COS ( XL*KR3 ) +V 1*S I N( XL*KR3 ) ) *XK3** 2 
*EXP( -XL*KI4 ) *( COS( XL*KR4 ) +V1*SIN( XL*KR4 ) ) *XK4**2 
=EXP( -XL*KI 1)*( COS( XL*KR 1 ) +V 1*SIN( XL*KR1) )*XK1**3 
=EXP( -XL*KI 2 ) *( COS ( XL*KR2 ) +V 1*SIN( XL*KR2 ) )*XK2**3 
=EXP( -XL*KI 3 ) *( COS( XL*KR3 ) +V1*S INC XL*KR3 ) ) *XK3**3 
sEXP(-XL*KI4)*(C0S(XL*KR4)+Vl*SIN(XL*KR4))*XK4**3 
n=(  F4*M3*P2 - F4*M3*P 1 -F4*M2*P3+F4*M2*P 1+F4*M 1*P3 -F4*M 1*F2 
-F3*M4*P2+F3*M4*P1+F3*M2*P4-F3V'M2*P1-F3*M1*P4+F3’'-M1*P2 
+F2*M4*P3-F2*H4*P1-F2*M3*P4+F2*II3*P1+F2*M1'VP4-F2*M1*P3 
-F1*M4*P3+F i*M4*P2+F 1*M3*P4-F1*M3*F2-F1*M2*P4+F 1*M2*F3 ) 
=(  F4*M3*P2 -F4*H2*P3-F3*M4*P2+F3*M2*P4+F2*M4*P3 -F2*M3*P4 ) 
/DEN 

-( -F4*M3*Pl+F4*Ml*P3+F3*M4*Pl-F3*tll*P4-Fl*M4*P3+Fl*M3*P4) 
/DEN 

=(  F4*M2*P 1 -F4*M1*P2-F2*M4*P1+F2*M1*P4+F 1*M4*P2-F1*M2*P4 ) 
/DEN 

=( -F3*M2*P 1+F  3*M 1 * P 2  +F 2*M 3 * P 1 -F2*M1*P3 -F1*M3*P2+F 1*M2*P3 ) 
/DEN 

S-V1*E*XI*( A1*XK1**3+A2*XK2**3+A3*XK3**3+A4*XK4**3) 


1 


IMPEDANCES 

IMF=2*FF/(V1*W) 
IMPR=REAL( IMP) 
IMPI=AIHAG( IMP) 


WRITE ( 15,1000)  F, IMPR, IMP I 
F=F+20.  0 
100  CONTINUE 

1000  FORMAT! IX, D9. 3,2X,D12.  5 ,2X,D12.  5) 

CLOSE! UN IT= 15) 

STOP 

END 

3.  THE  IMPEDANCE  OF  THE  FINITE  VISCOELASTIC  BEAM  USING 


TIMOSHENKO  BEAM  THEORY 


’  Vc  Vc  Vc  Vc  Vc  Vc  Vc  *  VcVc  Vc  Vc  Vc  V V  V?  Vc  ■ 


•  Vc  Vc  Vc  Vc  Vc  Vc  Vc  Vc  Vc  Vc  Vc  • 


THE  THREE  REDUCE  PROGRAMS  WERE  USED  FOR  EVALUATING  THE  IMPEDANCE 
OF  THE  FINITE  VISCOELASTIC  BEAM  USING  TIMOSHENKO  BEAM  THEORY. 


C 

C  THE  FIRST  REDUCE  PROGRAM  DEFINED  SHEAR  FORCE  AND  FOUR  BOUNDARY 
C  CONDITIONS. 

C 

YY:  =Y1*EXP( I*K1*X)+Y2*EXP( I*K2*X)+Y3*EXP( I*K3*X)+Y4*EXP( I*K4*X); 

PP1: =R1*Y1*£XP( I * K 1 *X ) +R 2 * Y 2 * E X P ( I*K2*X); 

PP2:  =R3*Y3*F.XP(  I*K3*X)+R4*Y4*EXP(  t*K4*X); 

PP: =PP1+PP2; 

MM: =DF( PP ,X); 

SS: =-(DF(YY,X)-PP); 

Y0:=SUB(X=0,YY)-Y; 

PO: =SUB( X=0 ,PP) -P; 

ML:  =SUB(X=L,NM); 

SL:  =SUB(.X=L,SS); 

BYE; 

C 

G  THE  SECOND  REDUCE  PROGRAM  SOLVED  CONSTANT  VARIABLES  (Yl,  Y2 , 

C  Y3  AND  Y4) . 

C 

YO:  =Y1+Y2+Y3+Y4-Y; 

P0:=R1*Y1+R2*Y2+R3*Y3+R4*Y4-P; 

ML: =M 1 * Y 1 +M  2 - Y 2 +M 3 * Y3+M4* Y4 ; 

SL:  =  S 1 *  Y 1 + S  2* Y2+ S  3* Y 3+S 4* Y4 ; 

SOLVE! LST(YO,PO, ML, SL),Y1,Y2,Y3,Y4); 

BYE; 

C 

C  THE  THIRD  REDUCE  PROGRAM  EVALUATED  THE  SHEAR  FORCE. 

P 

VJ 

YY:  =Y1*EXP( I*K1*X)+Y2*EXP( I*K2*X)+Y3*EXP( I*K3*X)+Y4*EXP( I*K4*X); 

PP1:  =R1*Y1*EXP( I * K 1 * X ) +R 2 * Y 2 * E X P ( I*K2*X); 

PP2:  =R3'VY3*EXP(  I*K3*X)+R4*Y4*EXP(  I*K4*X); 

PP: =PP1+PP2; 

S  S  S :  =  ■ -  KA*GG* A A* ( DF ( Y Y , X ) - P  P ) ; 

SS:  =SUB( X=0 , SSS) ; 

BYE; 


* 


THIS  FORTRAN  PROGRAM  DEVELOPED  THE  IMPEDANCE  OF  THE  FINITE 


VISCOELASTIC  BEAM  USING  TIMOSHENKO  BEAM  THEORY. 


DEFINED  VARIABLES; 


E  ; 

F  : 

G  : 

GRAV  : 
IMP  : 
KAR  : 
POISS  : 
ROH  : 
S3  : 

XB,  XH , XL 
XI  : 

XK1 ,  XK2 , 
XNETA 
VI 


COMPLEX  YOUNG'S  MODULUS  OF  THE  VISCOELASTIC  BEAM 
FREQUENCY  U HZ') 

COMPLEX  SHEAR  MODULUS  OF  THE  VISCOELASTIC  BEAM 

GRAVITY  (  32-6  IN/SEC 

IMPEDANCE  OF  THE  VISCOELASTIC  BEAM 

SHAPE  OF  CROSS  SECTION 

FOISSON’S  RATIO 

DENSITY  OF  THE  VISCOELASTIC  BEAM 
SHEAR  FORCE  (LBF) 

:  DIMENSIONS  OF  THE  BEAM  (WIDTH,  HEIGHT,  LENGTH) 
AREA  MOMENT  OF  INERTIA 
XK3 ,  XK4  :  WAVE  NUMBERS 

ENERGY  LOSS  FACTOR  OF  THE  VISCOELASTIC  BEAM 
IMAGINARY  (0,1) 


:‘V  rVc’jWrVc’VrjV’ 


Vr  Vr  “V  iY  '.V  Vr  Vc  Vc  Vc  *  Vr  V?  V?  ->V  rr  i  Y  3  V  -jY  ir  Vc  Vc  Vc  V : 


•  V?  “V  *  “Wr  VcVcVr  ■ 


COMPLEX  V 1 , G , E , XKN 1 , XKN2 , XKD , XK 1 , XK2 , XK3 , XK4 , R1 ,R2 ,R3 ,R4 
COMPLEX  Ml , M2 , M3 ,M4 ,S 1 ,S2 ,S3 , S4 ,DEN , A1 , A2 , A3 , A4 , SS , IMP 
REAL  F , KAR , PO I SS , GRA V , P I , XB , XL , XH , XG , ROH , XI , W , XLOG , XNETA 
REAL  KR1 , KR2 , KR3 , KR4 , KI 1 , KI 2 , KI 3 , KI4 , IMPR , IMPI 

0PEN(UNIT=15,FILE='TIM' ) 

CONSTANTS 

F=10. 0 
KAR=0. 83 
P0ISS=0. 45 
GRAV=386.  0 
PI=3. 1415926 
V1=CMPLX(0.  0,1.  0) 

DIMENSIONS 

XB=1  0 
XL=8. 0 
XH=0 . 35 


PROFERITIES 

RCH=0. 05505/GRAV 
X I = ( XB  •''XH 3 )  /  1 2 .  0 

FOUR  DIFFERENT  WAVE  NUMBERS 

DO  100  1=1,100 
V=2-VP  I*F 

XLOG=ABS( L0G10( F/62) ) 

XNETA=0. 65*EXP( -0. 52732*XLOG**l. 956) 
XNETA=0.  2 
XNETA=0.  5 


non  non 


XG=0. 0  0  0 0  2 5  0  3*F**  3 - 0 . 1752*F**2+457. 5883*F+29280 
C  XG=200000 

G=XG*CMPLX(  1.  O.XNETA) 

E=2*G*( 1+FOISS) 

XKN1=W**2*( KAR*G+E  1  /ROH 
XKN2=  SQRT(W**4*(  (  KAR*G+E)/R0H)**2 
&  -4*KAR*G*E*(V**4-KAR*G*XB*XH*W**2/ ( ROH*XI ) ) 

&  /  (  ROH**  2 ) ) 

XKD=2*KAR*G*E/( ROH**2 ) 

XK1=  CSQRT( (XKN1+XKN2)/XKD) 

XK2=-CSQRT( (XKN1+XKN2) /XKD) 

XK3=  CSQRT( fXKNl-XKN2) /XKD) 

XK4=-CSQRT((XKN1-XKN2)/XKD) 

KR1=REAL(XK1) 

KR2=REAL(XK2) 

KR3=REAL( XK3 ) 

KR4=REAL(XK4) 

K 1 1 = A I M  AG  ( XK 1 ) 

KI2=AIHAG(XK2) 

KI3=AIMAG(XK3) 

K I 4= A I MAG ( XK4 ) 

DEFINED  CONSTANTS  FROM  THE  THREE  REDUCE  PROGRAMS 

R 1=( ROH*W** 2 - KAR*G*XK 1**2 ) / ( V 1*KAR*G*XK 1 ) 

R2=( ROH* W** 2 - KAR*G*XK 2 **  2 ) / ( V1*KAR*G*XK2 ) 

R3=( ROH*W**2 -KAR*G*XK3**2 ) / ( V1*KAR*G*XK3 ) 

R4=( ROH*W**  2 - KAR*G*XK4**  2 ) / ( V 1*KAR*G*XK4 ) 

M 1=EXP ( -XL*KI 1 ) * ( COS ( XL*KR 1 ) +V 1*S IN( XL*KR 1 ) ) *V1 *XK 1*R 1 
M2=EXP( -XL*KI2)*( COS( XL*KR2 )+Vl*SIN( XL*KR2 ) )*V1*XK2*R2 
M3=EXP( -XL*KI 3 )*( COS( XL*KR3 ) +V 1*S IN( XL*KR3 ) ) *V1*XK3*R3 
M4=EXP( -XL*KI4) *( COS ( XL*KR4 ) +V 1*S IN( XL*KR4 ) ) *V1*XK4*R4 
S 1=EXP( -XL*KI 1 )*( COS( XL*KR1 )+Vl*SIN( XL*KR1) )*( R1 - V1*XK1 ) 
S2=EXP( -XL*KI2)*(COS(XD*KR2)+Vl*SIN( XL*KR2) )*(R2-V1*XK2 ) 

S  3=EXP ( -XL*KI 3 ) * ( COS ( XL*KR3 ) +V 1*S IN ( XL*KR3 ) ) *( R3 - V 1*XK3 ) 
S4=EXP( -XL*KI4 ) *( COS ( XL*KR4)+V 1*SIN( XL*KR4 ) ) *( R4 - V 1*XK4 ) 

DEN=(  S4*M3*R2-S4*M3*R1-S4*M2*R3+S4*M2*R1+S4*M1*R3-S4*M1*R2 
&  - S  3*M4*R2+S 3*M4*R 1+S3*M2*R4 - S3*M2*R 1 - S  3*M 1*R4+S 3*M 1*R2 

&  +S 2*M4*R3 - S  2*M4*R 1 - S  2*M3*R4+S 2*M3*R 1+S  2*M 1*R4 - S  2*M 1*R  3 

&  -S 1*M4*R3+S 1*M4*R2+S 1*M3*R4 -S 1*M3*R2 -S 1*M2*R4+S 1*M2*R3 ) 

A 1=(  S4*M3*R2 - S4*M2*R3 - S  3*M4*R2+S 3*M2*R4+S 2*M4*R3 - S 2*M3*R4 ) 

&  /DEN 

A2=( - S4*M3*R 1+S4*M 1 *R3+S3*M4*R 1 - S  3*M 1*R4 - S 1*M4*R3+S 1*M3*R4) 

&  /DEN 

A3=(  S4*M2*R 1 -S4*M 1*R2 - S2*M4*R 1+S 2*M 1*R4+S 1*M4*R2 -S 1*M2*R4 ) 

&  /DEN 

A4=( - S  3*M2*R 1+S3*M 1*R2+S2*M3*R 1 - S  2*M 1*R3 - S 1*M3*R2+S 1*M2*R3 ) 
h  /DEN 

SS=KAR*G*XB*XH*(  A 1*( R 1 -V1*XK1 )+A2*( R2-V1*XK2 )+A3*( R3 -V1*XK3 ) 
&  +A4*(R4-V1*XK4)  ) 

IMPEDANCES 

IMP=2*SS/(V1*W) 

IMPR=REALf IMP) 

IMPI=AIMAG( IMP) 
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WRITE! 15 , 1000)  F , IMPR , IMPI 
F=F+20. 0 
100  CONTINUE 

1000  FORMAT! IX, D9.  3,2X,D12. 5 ,2X,D12. 5) 

CLOSE(  UN'IT=  15) 

STOP 

END 

4.  THE  IMPEDANCE  OF  THE  FINITE  CONSTRAINED  LAVER  BEAM 
l ’SING  THE  SIXTH  ORDER  BEAM  THEORY 


THE  THREE  REDUCE  PROGRAMS  WERE  USED  FOR  EVALUATING  THE  IMPEDANCE 
OF  THE  INFINITE  VISCOELASTIC  BEAM  USING  BERNOULLI -EULER  BEAM 
THEORY. 


C 

C  THE  FIRST  REDUCE  PROGRAM  DEFINED  TRANSVERSE  FORCE  AND  TWO  BOUNDARY 
C  CONDITIONS. 

C 

YY1:  =Y1*EXP( I*K1*X)+Y2*EXP( I*K2*X); 

YY2:  =Y3*EXP( I*K3*X )+Y4*EXP( I*K4*X); 

YY3:  =Y5*EXP( I*K5*X)+Y6*EXP( I*K6*X); 

UU1:  =R1*Y1*EXP( I*K1*X)+R2*Y2*EXP( I*K2*X); 

UU2: =R3*Y3*EXP( I*K3*X)+R4*Y4*EXP( I*K4*X) ; 

UU3:  =R5*Y5*EXP( I*K5*X)+R6*Y6*EXP( I*K6*X) ; 

YYY: =YY1+YY2+YY3; 

UUU: =UU1+UU2+UU3; 

YO: =SUB(X=0 , YYY) -Y; 

PO:  =SUB(X=0 ,DF( YYY ,X) ) -P; 

UO: =SUB(X=0 ,UUU) -U; 

SL:  =SUB(X=L,(P P*DF( YYY , X . 3 ) +QQ*DF( YYY , X ) +RR*UUU ) ) ; 

ML:  =SUB(X=L , ( PP*DF( YYY , X , 2 ) +SS*DF( UUU , X) ) ) ; 

NL:  =SUB ( X=L , ( TT*DF( UUU , X) ) ) ; 

BYE; 

C 

C  THE  SECOND  REDUCE  PROGRAM  SOLVED  CONSTANT  VARIABLES  (Y4,  Y5 , 

C  AND  Y6) . 

C 

YO:  =Y1+Y2+Y3+Y4+Y5+Y6-Y; 

PO: =P1*Y1+P2*Y2+P3*Y3+P4*Y4+P5*Y5+P6*Y6-P; 

UO:  =R  1  *  Y 1 +R  2*  Y  2 +R  3*Y  3+R4* Y4+R5*Y  5 +R6* Y6 - U; 

SL: =Sl*Yl+S2*Y2+S3*Y3+S4*Y4+S5*Y5+S6*Yo; 

ML:  =M1*Y1+M2*Y2+M3*Y3+M4*Y4+M5*Y3+M6*Y6; 

NL:  =  N 1  *  Y 1 +N  2  *  Y  2 +N  3  Y  3 + N  4  *  Y  4 + N  5  *  Y  5 + N  6  Y  6 ; 

SOLVE! LST! SL , ML , NL) , Y4 , Y5 , Y6 ) ; 

BYE; 

C 

C  THE  THIRD  REDUCE  PROGRAM  SOLVED  VARIABLES  (Yl,  Y2  AND  Y3)  USING 

C  Y4,  Y5  AND  Y6.  THE  THIRD  REDUCE  PROGRAM  EVALUATED  THE  SHEAR  FORCE. 

C 

Y4:  =B1*Y1+C1*Y2+D1*Y3; 

Y3 : =B  2* Y 1+C  2* Y2+D2*Y 3 ; 

Y6:  =B3*Y1+C3*Y2+D3*Y3; 

YO: =Y1+Y2+Y3+Y4+Y5+Y6-Y; 


PO: =P1*Y1+P2*Y2+P3*Y3+F4*Y4+P5*Y5+P6*Y6-P; 
UO:  *R1*Y1+R2*Y2+R3*Y3+R4*Y4+R5*Y5+R6*Y6-U; 
SOLVE ( LST( YO , PO , UO ) , Y1 , Y2 , Y3 ) ; 

YYY1:  =Y1*EXF(  I*K1*X)+Y2*EXP(  I*K2*X); 

YYY2:  =Y3*EXP( I*K3*X) FY4*EXP( 1*K4*X); 

YYY3:  =Y5*EXP(  I*K5*X)+Y6*EXP(  I*K6*X); 

YYYY:  =YYY 1 +YYY2+YYY3; 

UUU 1 : =R1*Y1*EXP( I*K1*X)+R2*Y2*EXP( I*K2*X); 
UUU2: =R3*Y3*EXP( I*K3*X)+R4*Y4*EXP( I*K4*X) ; 
UUU3: =R5*Y5*EXP( I*K5*X)+R6*Y6*EXP( I*K6*X); 
UUUU:  =UUU 1+UUU2+UUU3 ; 

3  S  S :  =P P*DF ( YYYY ,  X ,  3 ) +QQ*DF( YYYY , X ) +RR*UUUU; 
SS:  =SCB(X=0  , SSS); 

BYE; 


THIS  FORTRAN  PROGRAM  DEVELOPED  THE  IMPEDANCE  OF  THE  CONSTRAINED 
LAYER  BEAM  USING  THE  6TH  ORDER  BEAM  THEORY. 


DEFINED  VARIABLES; 

El,  E3  :  COMPLEX  YOUNG'S  MODULUS  OF  THE  CONSTRAINED  LAYER 

F  :  FREQUENCY  (HZ) 

G  :  COMPLEX  SHEAR  MODULUS  OF  THE  VISCOELASTIC 

GRAV  :  GRAVITY  (386  IN/SEC**2) 

HI,  H3  :  HEIGHT  OF  THE  CONSTRAINED  LAYER 

H2  :  HEIGHT  OF  THE  VISCOELASTIC 

II,  12  :  AREA  MOMENT  OF  INERTIA 

IMP  :  IMPEDANCE  OF  THE  CONSTRAINED  LAYER  BEAM 

KAR  :  SHAPE  OF  CROSS  SECTION 

POISS  :  POISSON’S  RATIO 

ROH1  :  DENSITY  OF  THE  CONSTRAINED  LAYER 

ROH2  :  DENSITY  OF  THE  VISCOELASTIC 

SS  :  SHEAR  FORCE  ( LBF) 

XB,  XL  :  DIMENSIONS  OF  THE  BEAM  (WIDTH,  LENGTH) 

XK1 ,  XK2 ,  XK3,  XK4,  XK5 ,  XK6  :  WAVE  NUMBERS 

XNETA  :  ENERGY  LOSS  FACTOR  OF  THE  VISCOELASTIC  BEAM 

VI  ;  IMAGINARY  (0,1) 


■jV  *  V  V?  Vr  i:  Vr  ■>  c  *V  Vr  V?  Vr  Vr  *,V  Vr  Vc  “iV  -,V  V?  Vr  ?V  VrVr  Ve  Vc  Vr  ‘V  *V  Vc  ■> V  Vr  *5  r  Vr  Vc  Vf  Vr  Vr  “V  *  V  •jV  Vr  Vr  Vc  * V  Vr  Vr  Vr  Vc  iV  Vc  "jV  i V  'V  *,V  -.’r  Vc  Vr  Vr  -.V  Vr  /V  Vr  V  *.V  Vr  “V 


c 

COMPLEX*32  V 1 , G , BB , P 1 , P2 , P3 , P4 , P5 , F6 , R 1 , R2 , R3 , R4 , R5 , R6 , QQ , RR 
COMP  LEX*  3  2  S 1 ,  S2 ,  S3 ,  S4 ,  S5  ,  S6 ,  M 1 ,  M2 ,  M3  ,  M4  ,  M5 ,  M6  ,  N 1 ,  N2  ,  N3  ,  N'4 ,  N5  ,  N(3 
COMPLEX*32  B1,C1,D1,B2,C2,D2,B3,C3,D3, DENI , DEN2 , A1 , A2 , A3 
COMPLEX  B , D , XK 1 , XK2 , XK3 , XK4 , XE5 , XK6 , XXK 1 , XXK2 , XXK3 , IMP 

REAL* 16  F , P I , GRAV , XB , XL , XG , HI , H2 , H3 , E 1 , E3 , I 1 , 13 , DD , ROH 1 , ROH2 
REAL* 16  CC,XM,TL(30) ,T2(30) ,T3(30) ,01( 30 ) ,Q2( 30) ,Q3(3G) 

REAL*  1 6  XLOG ,  W ,  Z  ,  C ,  BR ,  B I ,  CR , DR ,  DI , F T 1 ,  FQ 1 ,  DTT 1 ,  DTQ 1 .  DQQ  l ,  DQT .1 
REAL* 16  FT 2 , FQ2 , DTT 2 , DTQ2 , D0Q2 , D0T2 , FT 3 , FQ3 , DTT 3 , DTQ 3 , DQQ3 , DQT3 
REAL* 16  XKR 1 , XKR2 , XER3 , XKI 1 , XKI2 , XKI3 , KR1 , KR2 , KR3 , KR4 , KR5 , XR6 
REAL* 16  KI1,KI2,KI3,KI4,KI5,KI6,PP,S3,TT, IMPR , IMP I 
REAL  XNETA 
C 

0PEN(UNIT=15,FILE='SAN' ) 

C 


I-B195  US  ANALYTICAL  AND  EXPERIMENTAL  STUDIES  OF  BEAM  WAVEGUIDE 
ABSORBERS  FOR  STRUCTURAL  DAHPING(U)  NAVAL  POSTGRADUATE 
SCHOOL  MONTEREY  CA  G  G  LEE  MAR  88 


UNCLASSIFIED 


F/G  13/18  1  NL 


CONSTANTS 


F=2010 

PI=3. 1415926 
GRAV=336.  0 
V1=CMFLX( 0.  0, 1-  0) 

DIMENSIONS 


XB=1.  0 
XL=8. 0 
H1=0. 0625 
H2=0.  065 
H3=0. 0625 

PROPERITIES 

El=10000000.  0 
E 3=10000000.  0 
I 1=XB*H1**3/ 12. 
I3=XB*H3**3/12. 
DD=E1*I1+E3*I3 
ROH1=0. 09832/GRAV 
R0H2=0. 03663/GRAV 
CC=Hl/2. +H2+H3/2. 

XM=XB* ( ( H 1+H3 ) *ROH 1+H2*R0H2 ) 


INITIAL  GUESS  WAVE  NUMBERS 

T 1 ( l)=-4.  31911362 
T2( l)=-0. 05139603 
T3( 1)=  3. 77697017 
Q 1 ( l)  =  -0.  33835599 
Q2( l)=-0. 05104052 
Q3(l)=-0.  27669665 

FIND  SIX  DIFFERENT  WAVE  NUMBERS  USING  NEWTON'S  METHOD 
SIX  DIFFERENT  WAVE  NUMBERS  ARE  XK1 ,XK2 ,XK3 ,XK4 ,XK5 ,XK6 

DO  100  1=1,100 
W=2*PI*F 

XLOG=ABS( LOG10(F/800) ) 

XNETA=1.  1*EXP( -0.  3906*XLOG**1. 53) 

XNETA=0.  5 
XNETA=1.  0 

XG=-0.  0000 125 8*F**2+0.  2472*F+74. 988 
XG=1000 

G=XG*CMPLX( 1.  0 ,XNETA) 

B  B=G*( E 1*H1+E3*H3 ) / ( H2*E 1*H1*E3*H3) 

Z=CC**2*E 1*H 1*E3*H3/ ( DD*(E 1*H1+E3*H3 ) ) 

B=BB*(  1.  0+Z) 

C=-XM*W**2/DD 
D=-BB*XM*W**2/DD 
BR=REAL( B) 

BI=AIMAG( B) 

CR=C 
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DR=REAL(Dt 
DI=AIMAGf D) 

DO  200  J=1 , 15 

FT  1=T  1  ( J  'l  **3+BR*T  1  (  J  V **2+(  -  3*Q  llJ)  **2  -  2*B  I  *Q  1(  J'i+CR)  *T  1(  J 1 
Sc  +( -BR*Q1(  J)**2+DR) 

FQ1=01(J)**3+BI*Q1( J)**2+(-3*Tl(J)**2-2*BR*Tl(J)-CR)*Ql( J) 
&  +(-BI*Tl(J)**2-DI) 

DTT1=3*T1( J)**2+2*BR*T1( J) 

&  +( - 3*0 1 ( J ) ** 2 - 2*  B I*Q 1 ( J ) +CR ) 

DTQ 1=( -6*T1 ( J) - 2*BR )*Q 1 ( J ) - 2*B I*T1 ( J) 

DQQ 1=3*Q 1 ( J ) **  2+2* B I*Q 1 ( J ) 

Sc  +(  -3*T1(  J)**2 -2*BR*T1(  J)  -CR3 

DOT 1=( - 6 *Q 1 ( J ) - 2* B I ) *T 1 ( J ) - 2*BR*0 1 ( J ) 

T 1 ( J+ 1 ) =T1 ( J) - (  DQQ 1*FT1 -DTQ 1*FQ l ) / ( DTT 1*DQQ 1 -DTQ 1*DQT 1 ) 

Q 1 ( J+  n  =Q 1 ( J ) - ( -DQT 1*FT1+DTT1*FQ 1 ) / ( DTT 1*D0Q 1 - DTO l*DOT 1 ) 
FT2=T2 ( J ) **3+BR*T2( J)**2+( - 3*Q2( J ) **2 - 2*B I*Q2 ( J ) +CR ) *T2 ( J ) 
&  +( - BR*Q2( J )**2+DR ) 

F  Q  2 =Q  2  (  J )  **  3 + B I  *  Q  2  (  J )  **  2 + ( -3*T2( J)**2-2*BR*T2( J) -CR)*Q2( J) 
&  +(-BI*T2( J)**2-DI) 

DTT2=3*T2( J)**2+2*BR*T2( J) 

Sc  +(  -  3 * Q 2  (  J  )  ** 2  -  2 * B I  * Q 2  f  J  )  +C R ) 

DTQ2=( - 6*T2 ( J ) - 2*B  R ) *Q2 ( J ) - 2*B I*T2 ( J ) 

DQQ2=3*Q2( J)**2+2*BI*Q2( J) 

Sc  +(  -3*T2(  J)**2-2*BR*T2(  J)  -CR) 

DQT2=( -6*Q2( J) -2*B I )*T2( J) -2*BR*Q2C  J) 

T2( J+l )=T2( J) - (  DQQ2*FT2 -DTQ2*FQ2) / ( DTT2*DQQ2-DTQ2*DOT2 ) 
Q2( J+l )=Q2( J) - ( -DQT2*FT2+DTT2*FQ2 ) / ( DTT2*DQQ2 -DTQ2*DQT2 ) 
FT3=T3(  J)**3+BR*T3(  J)**2+(  -3*03(  J)**2-2*BI*Q3(  J)+CR)*T3('  J) 
&  +( -BR*Q3( J)**2+DR) 

FQ3=Q3( J)**3+BI*Q3f J)**2+( -3*T3( J)**2-2*BR*T3( J)-CR)*Q3( J) 
&  +(-BI*T3(J)**2-DI) 

DTT3=3*T3 ( J ) **2+2*B  R*T3 ( J ) 

Sc  +(  -  3 *Q 3 ( J )  ** 2  -  2* B I *Q 3  ( J )  +CR ) 

DTQ3=( - 6*T3 ( J ) - 2*BR ) *Q3 ( J ) - 2*B I*T3 ( J ) 

DQQ3=3*Q3(  J)**2+2*BI*Q3(.  J) 

Sc  +(  -3*T3(  J)**2-2*BR*T3(  J)  -CR) 

DQT3=( -6*Q3( J) -2*BI)*T3( J) -2*BR*Q3( J) 

T3 ( J+ 1 )=T3( J ) - (  DQQ 3* FT 3- DTQ3*FQ3 ) / ( DTT3*DQQ3 -DTQ 3*DQT3 ) 
Q3( J+l )=03( J) -( -DQT3*FT3+DTT3*FQ3 ) / ( DTT3*DQQ3 -DTQ3*DQT3 ) 
200  CONTINUE 

XKR1=  Tl( J) 

XKR2=  T2(J) 

XI<R3=  T3(  J) 

XKI 1=  Q 1 ( J ) 

XKI2=  Q2( J) 

XKI 3=  Q3( J) 

XXK  l=c:  1PLX(  XKR 1 ,  XKI  1 ) 

XXK2=CIIPLXfXKR2,XKI2) 

XXK3=CMPLX( XKR 3 , XKI3 ) 

XKI®  CSQRT(XXKl) 

XK2=-CSQRT( XXK 1 ) 

XK3=  CSQRT(XXK2) 

XK4=-CSQRT(XXK2) 

XK3®  CSQRTf  XXK 3) 

XK6=-CSQRT( XXK3) 

KR1=  REAL(XKl) 


KR2=  REAL(XK2) 

KR3=  REAL(XK3) 

KR4=  REAL(XK4) 

KR5=  REAL(XK5) 

KR6=  REAL(XK6) 

KI1=AIMAG(XK1) 

K 1 2=A I MAG( XK  2 ) 

KI 3=AIMAG( XK3 ) 

K 1 4= A I M  AG  (  XK4  'i 
K 1 5 = A I MAG ( XK5 ) 

KI6=AIMAG(XKo) 

DEFINED  CONSTANTS  FROM  THE  THREE  REDUCE  PROGRAMS 

P1=V1*XK1 

P2=V1*XK2 

P3=V1*XK3 

P4=V1*XK4 

P5=V1*XK5 

P6=V1*XK6 

R 1=V1*BB*Z*DD*XK 1 / ( ( XK1**2+BB ) *E3*H3*CC ) 

R2=V  1*BB*Z*DD*XK2/ ( ( XK2**2+BB )*E3*H3*CC ) 
R3=Vl*BB*Z*DD*XK3/(  ( XK3**2+BB )*E3*H3*CC) 
R4=Vl*BB*Z*DD*XK4/((XK4**2-f-BB')*E3*H3*CC) 

R5=V 1*B  B*  Z*DD*XK5 / ( ( XK5**2+BB )*E3*H3*CC ) 
R6=V1*BB*Z*DD*XK6/  ( ( XK6**2+BB )*E3*H3*CC ) 

PP=DD 

QQ=-DD*BB*Z 

RR=DD*BB*Z*( E 1*H 1+E3*H3 ) / ( CC*E 1*H 1 ) 

SS=-CC*E3*H3 

TT=E3*H3 

S 1=( - V 1*PP*XK 1**3+V 1*QQ*XK 1+RR*R 1 ) *EXP( - XL*KI 1 ) 

&  *( COS(XL*KR 1 ) +V 1*S IN( XL*KR 1 ) ) 

S  2=( - V 1 *pp*xK2**3+V 1*QQ*XK2+RR*R2 ) *EXP ( -XL*K 1 2 ) 

&  *( COS( XL*KR2)+V1*SIN( XL*KR2) ) 

S3=( -V l*pp*XK3**3+V 1*QQ*XK3+RR*R3 ) *EXP( -XL*KI 3 ) 

6c  * (  COS ( XL*KR3  )  +V  1*S  I N(  XL*KR3  )  ) 

S4=(  - V  1*PP*XK4**  3+V  1*QQ*XK4+RR*R4 ) *EXP(  -XL*KI4 ) 

6c  *(  COS (  XL*KR4 )  +V  1*S INC  XL*KR4 ) ) 

S5=( -V1*PP*XK5**3+V1*QQ*XK5+RR*R5)*EXP( -XL*KI5 ) 

&  *(C0S(XL*KR5)+V1*SIN(XL*KR5) ) 

S  6=( - V 1*PP*XK6**3+V 1 *QQ*XK 6 +RR* R 6 ) *EXP ( -XL*K 1 6 ) 

&  *( COSC XL*KR6 ) +V 1*S IN( XL*KR 6 ) ) 

M 1=( V 1*SS*R1*XK 1 -PP*XK1**2 )*EXP( -XL*KI 1 ) 

6c  *(  COS(XL*KR  1 )  +V  1*SIN(  XL*KR  1 ) ) 

M2=( V1*SS*R2*XK2-PP*XK2**2 )*EXP( -XL*KI 2 ) 

6c  *(  COS ( XL*KR2 )  +V  1*S  IN(  XL* KR 2 ) ) 

M3=( V1*SS*R3*XK3 -PF*XK3**2 )*EXP( -XL*KI3 ) 

6c  *(  COS(  XL*KR3)+V1*SINC  XL*KR3I ) 

M4=( V1*SS*R4*XK4-PP*XK4**2)*EXP( -XL*KI4) 

6c  *(  COS(  XL*KR4)+V1*SIN( XL*KR4)  ) 

M 5  - ( V 1  *  S  S *  P.  5  * XK5  -  P  P *XK 5  **  2. )  * E  X F  (  -  X L*  K 1 5 ) 

5c  *(  COS(  XL*KR3 )  +V1*SIN(  XL*KR5 ) ) 

M6=(  V1*SS*R6*XK6 -PP*XK6**2 )*EXP( -XL*KI6 ) 

&  *( COS( XL*KR6 ) +V1*SIN( XL*KR6 ) ) 

N1=V1*TT*R1*XK1*EXP( -XL*KI 1) 


(n  <?■>  (r  fr  fr 


&  *( COS(  XL*KR1)  +V1*SIN(  XL*KR  1 )  'l 

N2=V 1*TT*R2*XK2*EXP( -XL*KI 2 ) 

&  *(C0S(XL*ER2)+V1*SIN(XL*KR2)) 

N3=V  1*TT*R3*XK3*EXF (  -  XL-KI 3  'i 
&  *f COS( XL*KR3 ) +V1*SIN( XL*KR3) ) 

S4=vi*TT*R4*xK4*EXP( -XL*KI4) 

&  *( COS ( XL*KR4 ) +V 1*S IN( XL*KR4) ) 

N 5 = V 1 *TT* R 5 * XK 5 * E X P ( -XL*KI5  > 

&  *( COS(XL*KR5  'I  +V1*SIN(  XL*KRS ) ) 

N6=V1*TT*R6*XK6*EXP( -XL*KI6 ) 

&  *( COS( XL*KR6)+V1*SIN( XL*KR6 ) ) 

B1=(-N6*M5*S1+N6*M1*S5+N5*M6*S1-N5*M1*S6-N1*M6*S5+N1*!!5*S6) 
&  /(N6*M5*S4-N6*M4*S5-N5*M6*S4+N5*M4*S6+N4*M6*S5-N4*M3*S6) 

Cl=( -N6*M5*S2+N6*M2*S5+N5*M6*S2-M5*M2*S6-N2*M6*S5+N2*M5*S6) 
&  /(N6*M5*S4-N6*M4*S5-N5*M6*S4+N5*M4*S6+N4*M6*S3-N4*M5*S6) 

Dl=f -N6*M5*S3+N6*M3*S5+N5*M6*S3-N5*M3*S6-N3*M6*S5+N3*M5*S6 > 
&  /(N6*M5*S4-N6*M4*S5-N5*M6*S4+N5*M4*S6+N4*M6*S5-N4*M5lVS6) 

B2=(N6*M4*SI-N6*M1*S4-N4*M6’>S1+N4*M1*S6+N1*M6*S4-N1*M4*S6 ) 

&  /(N6*M5*S4-N6*M4*S5-N5*M6*S4+N5*M4*S6+N4*M6*S3-N4*!J5iVS6) 

C2=(N6*M4*S2-N6*M2*S4-N4*M6*S2+N4*M2*S6+N2*M6*S4-N2*M4*S6) 

&  / ( N6*M5*S4 -N6*M4*S5 -N5*M6*S4+N5*M4*S6+N4*M6*S5 -N4*M5*S6 ) 

D2=(N6*M4*S3-M6*M3*S4-N4*M6*S3+N4*M3*S6+N3*M6*S4-N3*M4*S6) 

&  / ( N6*M5*S4 -N6*M4*S5 “N5*M6*S4+N5*M4*S6+N4*M6*S5 -N4*f  15*S6 ) 

B3=( -N5*M4*S 1+N5*M 1*S4+N4*MS*S 1 -N4*M 1*S  5 -N 1*M5*S4+N 1**14* S3) 
Sc  / ( N6*M5*S4-N6*M4*S5 -N5*M6*S4+N5*M4*S6+N4*H6*S5 -N4*M5*S6 ) 

C3=( -N5*M4*S2+N5*M2*S4+N4*M5*S2-N4*M2*S5 -N2*M5*S4+N2*M 4*S5 ) 
&  /  ( N6*M5*S4-N6*M4*S5  -N5*M6*S4+N5*M4*S6+N4*M6*S5  -N4*M5*S6 ) 

D3=(-N5*M4*S3+N5*M3*S4+N4*M5*S3-N4*M3*S5-N3*M5*S4+N3*M4*S5) 
&  /(N6*M5*S4-N6*M4*S5 -N5*M6*S4+N5*M4*S6+N4*M6*S5 -N4*M5*S6) 

DEN1=( R6*P5 -R5*P6)*( D3*C2*B 1+D3*C2 -D3*B2*C1-D3*B2-C3*D2*B 1 
-C3*D2+C3*B2*D1+C3*B2+B3*D2*C1+B3*D2 
*  B  3*C  2*  D 1 - B 3*C  2 ) 

+(R4*P6-R6*P4)*(D3*C2*B1-D3*B2*C1-D3*C1+D3*B1-C3*D2*B1 
+C3*B2*D  1+C  3*D  1  -C3*B  1+B  3*D2*C  1  -  B  3*C  2*D  1 
-B3*D1+B3*C1) 

+( R5*P4-R4*P5 ) *( D3*C2*B 1 -D3*B2*C 1 -C3*D2*B 1+C3*B2*D 1 
+B3*D2*C  1  -B3*C2*D  ]  +D2--C 1  -D2*B  1  -C2*D 1 
+C2*B 1+B2*D 1 -B2*C 1 ) 
DEN2=(R3*P6-R6*P3)*(C3*B2+C3*B1+C3-B3*C2-B3*C1-B3) 

&  +(R6*P2-R2*P6)*(D3*B2+D3*B1+I)3-B3*D2-B3*D1-B3) 

Sc  +(R1*P6-R6*P1)*(D3*C2+D3*C1+D3-C3*D2-C3*D1-C3) 

a  +(R5*P3-R3*P5)*(C3*B2-B3*C2-C2*B1-C2+B2*C1+B2) 

Sc  +(R2*P5-R5*P2)*(D3*B2-B3*D2-D2*B1-D2+B2*D1+B2) 

5c  +( R5*P  1-R1*P5  )*(  D3*C2  -C3*D2  -D2"C  1 -D2+C2--D1+C2 ) 

&  +(R4*P3-R3*P4)*(C3*B1-B3*C1+C2*B1-B2*C1-C1+B1) 

6-  +( R  2  *  P  4  -  R  4  *  P  ^ )  *  ( D3---B  1-B3*D1+D2*B1 -B2*D1 -Dl+Bl) 

5c  +( R4*P  1  -R1*P4  )*(  D3*C  1  -C3*D  1  +D2*C  1  -C2*D  1  -D  1+C  1 ) 

5c  +( R3*P2-R2*P3)*(  B3+B2+B 1+1 ) 

5c  +(R1*P3-R3*P1)*(C3+C2+C1  +  1) 

5c  +( R2'VP1-R1"  P2 )'”( D3+D2+D 1  +  1 ) 

A 1=(  R6*(  P5*D3*C2 - P5*C 3*D2+P4*D3*C 1 -P4*C3*D 1 - P3*C3+P2*D 3 ) 

5c  +  R5*(  -P6*D3*C2+P6*C3*D2+P4*D2*C1-P4*C2*D1-P3--C2+P2*D2) 

5c  +  R4*( -P6*D3*C1+P6*C3*D1-P5*D2*C1+P5*C2*D1-P3*C1+P2*D1) 

5c  +  R3''r(  P6*C3+P5*C2+P4*C1+P2) 

5c  +  R2*(-P6*D3-P3*D2-P4*D1-P3)  ) / ( DEN1+DEN2 ) 

A2=(  R6*( -P5*D3*B2+P5*B3*D2-P4*D3*B1+P4*B3*D1+P3*B3-P1*D3) 


o  o  o 


St  +  R5*(  P6*D3*B2-P6*B3*D2-P4*D2*B1+P4*B2*D1+P3*B2-P1*D2) 

St  +  R4*(  P6*D3*B1-P6*B3*D1+P5*D2*B1-F5*B2*D1+P3*B1-P1*D1) 

St  +  R  3  *  (  -  P  6  *  B  3  -  P  5  *  B  2  -  F  4  *  B 1  -  P 1 1 

St  +  Rl*(  P6*D3+P3*D2+P4*D  1+P3  )  ) /(  DEN1+DEN21 

A3=(  R6*(  P5*C3*B2-P5*B3*C2+P4*C3*Bi-P4*B3*Cl-P2*B3+Pl*C3) 

St  +  R5*( -P6*C3*B2+P6*B3*C2+P4*C2*B1-P4*B2*C1-P2*B2+P1*C2) 

St  +  R4*(-P6*C3*B1+P6*B3*C1-P5*C2*B1+F5*B2*C1-P2*B1+P1*C1) 

St  +  R2*(  P6*B3+P5*B2+P4*B1+P1) 

St  +  R 1* (  - P6*C 3 - P5*C 2 - P4*C  1  - P 2 )  )/(DENl+UrN2) 

IMPEDANCES 

IMP=2*(  ( R4*RR+V1*QQ*XK4-V1*PP*XK4**3)*(D1*A3+C1*A2+D  1*A1. ) 

St  +  ( R5*RR+V 1*QQ*XK3 -V 1*PP*XK5**3 )*f D2*A3+C2*A2+B2*A 1 ) 

St  +  ( R6*RR+V 1*QQ*XK6 - V 1*PP*XK6**3)*( D3*A3+C3*A2+B3*A1 ) 

&  +  RR*(R3*A3+R2*A2+R1*A1) 

S  +  V1*Q0*(XK3*A3+XK2*A2+XK1*A1) 

St  -V1*PP*(K3**3*A3+K2**3*A2+K1**3*A1)  )/(Vl*W) 

IMPR=REAL(  IMP 'i 
IMPI=AIMAG( IMP) 

WRITE ( 15 , 1000)F, IMPR, IMPI 
F=F-20 
Tl( 1)=T1( J) 

T2( 1 )=T2( J) 

T3( 1)=T3( J) 

QIC 1)=Q1( J) 

Q2( 1)=Q2( J) 

03 ( 1)=Q3( J) 

100  CONTINUE 

1000  FORMAT! IX, D9.  3,2X,D12.  5 ,2X,D12.  5) 

CLOSE( UNIT=15 ) 

STOP 

END 


IMPEDANCES  VS.  FREQUENCY 


APPENDIX  B. 


Figure  49. 


Real  part  of  the  20"  viscoelastic  beam  waveguide  absorber  impedances 
at  the  center  of  the  beam. 
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Figure  52.  Modal  damping  factors  vs.  frequencv  of  tiie  test  plate  with  a  1(> 
lastic  beam  waveguide  absorber  at  location  1  and  without. 
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Figure  53.  Tlie  driving  point  frequency  response  of  the  test  plate  with  a  20"  con¬ 


strained  layer  beam  waveguide  absorber  (dashed)  at  location  !  and 
without  ( solid). 


MODAL  DAMPING  VS.  FREQUENCY 
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Figure  55.  The  driving  point  frequency  response  of  the  test  plate  u it li  a  16"  viscoe¬ 
lastic  beam  waveguide  absorber  (dashed)  at  location  2  and  without 
(solid). 
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MODAL  DAMPING  VS.  FREQUENCY 


Figure  56.  Modal  damping  factors  vs.  frequency  of  the  test  plate  nitli  a  16"  viscoe 
lastic  lieani  \\a\eguide  absorber  at  location  2  and  without. 


Figure  57.  The  driving  point  frequency  response  of  the  test  plate  mtli  a  lt>'  con 
strained  layer  beam  waveguide  absorber  (dashed)  at  location  2  am 
without  (solid). 
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Figure  59.  The  driving  point  frequency  response  of  the  test  plate  with  a  20"  con¬ 
strained  layer  beam  waveguide  absorber  (dashed)  at  location  2  and 
without  (solid). 
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Figure  61.  The  driving  point  frequency  response  of  the  test  plate  with  a  16"  viscoe¬ 
lastic  beam  waveguide  absorber  (dasiied)  at  location  3  and  without 
(solid). 
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Figure  6<S.  Modal  damping  factors  vs.  frequency  of  the  test  plate  with  a  2lr  con 


strained  laver  beam  waveguide  absorber  at  location  3  and  without. 
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